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Abstract — This paper presents a coordination algorithm for 
mobile autonomous robots. Relying upon distributed sensing the 
robots achieve rendezvous, that is, they move to a common 
location. Each robot is a point mass moving in a nonconvex 
environment according to an omnidirectional kinematic model. 
Each robot is equipped with line-of-sight limited-range sensors, 
i.e., a robot can measure the relative position of any object 
(robots or environment boundary) if and only if the object is 
within a given distance and there are no obstacles in-between. 
The algorithm is designed using the notions of robust visibility, 
connectivity-preserving constraint sets, and proximity graphs. 
Simulations illustrate the theoretical results on the correctness 
of the proposed algorithm, and its performance in asynchronous 
setups and with sensor measurement and control errors. 

Index Terms — Multi-robot coordination. Cooperative control. 
Distributed algorithm. Mobile robot. Visibility, Nonlinear systems 
and control 



I. Introduction 

Multi-agent robotic systems have been receiving increasing 
attention in recent times. This is due in no small part to the 
remarkable advances made in recent years in the development 
of small, agile, relatively inexpensive sensor nodes with mobile 
and networking capabilities. These systems are ultimately 
intended for a wide variety of purposes such as search and 
rescue, exploration, environmental monitoring, location-aware 
computing, and the maintaining of structures. In this paper 
we design algorithms to steer a group of robots to a common 
location, or rendezvous, in a nonconvex environment. 

The rendezvous problem is a fundamental motion coordina- 
tion problem for collections of robots. In its essence, it is the 
most basic formation control problem and can be used as a 
building block for more sophisticated behaviors. It is related to 
the classic consensus problem in distributed algorithms. This 
problem is, here, tackled in a fully distributed manner, i.e., our 
robots do not have any global knowledge about the position 
of all other robots, do not have any global knowledge about 
the environment, and do not share a common reference frame. 
The information that is available to an individual robot is only 
what is provided by a local "visibility sensor." In other words, 
a robot can measure the relative position of a second robot if 
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and only if the robots are visible to each other in the nonconvex 
environment and lie within a given distance of each other. 

The literature on multirobot systems is very extensive. 
Examples include the survey in [1] and the special issue [2] 
of the IEEE Transaction on Robotics and Automation. Our 
multi-robot model is inspired by the literature on networks of 
mobile interacting robots: an early contribution is the model 
proposed in [3] consisting of a group of identical "distributed 
anonymous mobile robots" characterized as follows. Each 
robot completes repeatedly a cycle of operations: it senses 
the relative position of all other robots, elaborates this in- 
formation, and moves. In this early work, the robots share 
a common clock. A related model is presented in [4], where 
the robots evolve asynchronously, have limited visibility, and 
share a common reference frame. For these types of multirobot 
systems, the "multi-agent rendezvous" problem and the first 
"circumcenter algorithm" have been introduced in [5]. This 
circumcenter algorithm has been extended to various asyn- 
chronous strategies in [6], [4], where rendezvous is referred 
to as the "gathering" problem. The circumcenter algorithm has 
been extended beyond planar problems to arbitrary dimensions 
in [7], where its robustness properties are also characterized. 
None of these previous works considers the problem in non- 
convex environments with line-of-sight visibility sensors. 

We conclude the literature review by mentioning that forma- 
tion control and rendezvous problems have been widely inves- 
tigated with different assumptions on the inter-agent sensing. 
For example, a control law for groups with time-dependent 
sensing topology is proposed in [8]; this and similar works, 
however, depend upon a critical assumption of connectivity 
of the inter-agent sensing graph. This assumption is imposed 
without a model for when two robots can detect and measure 
each other's relative position. In this paper, we consider 
position-dependent graphs and, extending to visibility sensors 
a key idea in [5], we show how to constrain the robots' motion 
to maintain connectivity of the inter-agent sensing graph. 

Next, we describe the essential details of our model. We 
consider a group of robotic agents moving in a nonconvex 
environment without holes. We assume each robot is modeled 
as a point mass. We assume that each robot is equipped 
with an omnidirectional limited-range visibility sensor, the 
nomenclature is adopted from [9, Section 11.5]. Such a sensor 
is a device (or combination of devices) that determines within 
its line of sight and its sensing range the following quantities: 
(i) the relative position of other robots, and (ii) the relative 
position of the boundary of environment. By omnidirectional 
we mean that the field-of-vision for the sensor is 27r radians. 
Examples of visibility sensors are scanning laser range finders 
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Fig. 1. Execution of the Perimeter Minimizing Algorithm de- 
scribed in Section IV-AI on a group of robots distributed in a polygon, Q, 
shaped like a typical floor plan. The graph shown in the left-most figure is 
the r-range visibility graph Gr-vis,Q, (see Section llTl. 



with accurate distance measurements at high angular densitylj 
time-of-flight range cameras]! and optical depth sensors based 
on structured light systems, e.g., see [10]. The range data 
obtained from the sensors can be processed to obtain a 
geometric representation of the area visible from a robot, e.g., 
see [11]. We do not directly address in this work issues related 
to feature extraction from range data. We assume that the 
algorithm regulating the robots' motion is memoryless, i.e., 
we consider static feedback laws. Given this model, the goal 
is to design a discrete-time algorithm which ensures that the 
robots converge to a common location within the environment. 
See Figure [T| for a graphical description of a simulation. 

This paper's main contribution is a novel provably correct 
algorithm that achieves rendezvous in a nonconvex planar 
environment among robots with omnidirectional range-limited 
visibility sensors. Rendezvous is achieved among all robots if 
the inter-agent sensing graph is connected at the initial time 
or becomes connected at any time during the evolution. 

The technical approach contains a number of novel con- 
tributions and proceeds as follows. First, we review a few 
useful geometric notions [12], such as robust visibility [13] and 
proximity graphs [14], and introduce various novel visibility 
graphs. Second, to maintain connectivity during the system 
evolution, we design novel constraint sets that (i) ensure that 
the visibility between two robots is preserved, and (ii) are 
upper semicontinuous or closed maps of the robots' positions. 
Third, based on a discussion on visibility graphs, we define 
a new proximity graph, called the locally-cliqueless visibility 
graph, which contains fewer edges than the visibility graph, 
and has the same connected components. This construction is 
useful in the connectivity maintenance problem. Fourth, we 
provide a careful analysis of the algorithm we propose. As 
novel Lyapunov function, we consider the perimeter of the 
relative convex hull of the robot positions. The main theorem 
is proved via our recent version of the LaSalle Invariance 
Principle for set- valued maps [7]. Fifth and final, extensive 
simulations validate our theoretical results and establish the 
convergence of our algorithm under more realistic assumptions 
than the ones adopted in the theoretical analysis: our algorithm 
performance is still adequate assuming asynchronous agent 
operation, noise errors in sensing and control, or finite-size 
disk robots. 

'E.g., the Hokuyo URG-04LX, see |http: / / www . hokuyo-aut . jp| 
and the Sick S2000, see http: / / www . sick . com| 

^E.g., the SwissRanger SR-3000, see |http : / /www . swis granger . ch| 



The paper is organized as follows. Section HI] contains some 
useful geometric notions. In Section [III] we model robots 
with visibility sensors and we introduce the rendezvous and 
connectivity maintenance problems. In Section |IV] we intro- 
duce constraint sets and locally-cliqueless visibility graphs. 
In Section [V] we propose the Perimeter Minimizing 
Algorithm for the rendezvous problem. Numerical simula- 
tions are presented in Section |VT] Additional analysis results 
and all proofs are presented in Appendices ID HIIl 

II. Geometric notions 

In this section we introduce some useful geometric notions. 
We begin by reviewing some standard notation. Let Z>o, M, 
M>o, and M>o denote the sets of nonnegative integer, real, 
nonnegative real, and positive real numbers, respectively. For 
p G and r G M>o, let B{p, r) denote the closed ball 
centered at p of radius r. Given a bounded set X C M^, let 
co{X) denote the convex hull of X, and let CC{X) denote the 
circumcenter of X, i.e., the center of the smallest-radius circle 
enclosing X. For p, q G M^, let ]p, q[= {\p + (1 - A)(jf | < 
A < 1} and [p,q] = {\p + (1 - A)(j | < A < 1} denote 
the open and closed segment with extreme points p and q, 
respectively. Let \X\ denote the cardinality of a finite set X in 
M^. Given a set of points X C M^, and another point p G M^, 
let dist(p, X) denote the Euclidean distance of p to the set X. 
The diameter diam(X) of a compact set X is the maximum 
distance between any two points in X. 

Now, let us turn our attention to the type of environments 
we are interested in. Given any compact and connected subset 
Q of M?, let dQ denote its boundary. A point q of dQ is 
strictly concave if for all e > there exists qi and q2 in 
B{q, e) n dQ such that the open interval q2[ is outside Q. 
A strict concavity of dQ is either an isolated strictly concave 
point or a connected set of strictly concave points. According 
to this definition, a strict concavity is either an isolated point 
(e.g., points ri and r2 in Figure |2]i or an arc (e.g., arc ai in 
Figure |2]i. Also, any strictly concave point belongs to exactly 
one strict concavity. 

Definition II.l (Allowable environments) A set Q c is 

allowable if 

(i) Q is compact and simply connected; 

(ii) dQ is continuously dijferentiable except on a finite 
number of points; 

(iii) dQ has a finite number of strict concavities. 

Recall that, roughly speaking, a set is simply connected if it 
is connected and it contains no hole. A particular case of the 
environment described above is a polygonal environment, the 
concavities being the reflex vertice^ of the environment. 

At almost all strictly concave points v, one can define the 
tangent to dQ. (Here, the wording "almost all" points means 
all except for a finite number) At all such points v, the internal 
tangent half-plane Hq{v) is the half-plane whose boundary is 
tangent to dQ at v and whose interior does not contain any 
points of the concavity; see Figure |2] 

'a vertex of a polygon is reflex if its interior angle is strictly greater than tt. 
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Fig. 2. An allowable environment Q: the closed arc ai and the isolated 
points ri , r2 are strict concavities, v' is a point on ai where the slope of 
dQ is defined. Hq (d' ) is the half-plane with the tangent to dQ at v' as the 
boundary and the interior in the direction of the an'ow. v" is a point on ai 
where the slope of dQ is not well-defined. In this case, we define the tangent 
to be the one shown in the plot. Hq(v") is the half-plane with the tangent 
to dQ at v" as the boundary and the interior in the direction of the arrow. 



A point g g Q is visible from p G Q if [p, q\ C Q. The 
visibility set V{p) C Q from p G Q is the set of points 
in Q visible from p. This notion can be extended as follows 
(see [13]): 

Definition II.2 (Robust visibiUty) Take e > and Q cR"^. 

(i) The point q £ Q is e-robustly visible from the point 
peQ ifUg,^ip^^gjB{q',e) C Q. 

(ii) The e-robust visibility set V{p, e) C Q from p £ Q is 
the set of points in Q that are e-robustly visible from p. 

(iii) The e-contraction Qe of the set Q is the set {p G 
Q I \\p - q\\ > e for all q € dQ}. 

These notions are illustrated in Figure [3] Loosely speaking, 
two points p, q are mutually 0-robustly visible if and only if 
they are mutually visible. We present the following properties 
without proof in the interest of brevity. 




Fig. 3. Robust visibility notions. Q is the outer polygonal environment; the 
e-contraction Qf is the region with the curved boundary and containing the 
point p; the visibility set V(p) is the region shaded in light gray; the e-robust 
visibility set V(p, e) is the region shaded in darker gray. Note that the isolated 
concavities of Q give rise to strictly concave arcs in Qf . 



Lemma II.3 Given an allowable environment Q and e > 0, 
the following statements hold: 

(i) q G Q is e-robustly visible from p G Q if and only if 

\p, q] c Qe,- 

(ii) if e is sufficiently small, then is allowable; 



(iii) all strict concavities of dQ^ have non-zero length and 
are continuously differentiable. 

Remarks II.4 (i) In light of Lemma III.Bf ii). in what fol- 
lows we assume that e is small enough for to be 
connected and therefore allowable, 
(ii) Robust visibility is a useful concept in many practically 
meaningful ways. For example, according to this notion, 
points are visible only if they are at least at a distance 
e from the boundary. This is useful when an object is 
arbitrarily close to the boundary and is indistinguishable 
from the boundary itself. Additionally, the parameter e 
might be thought of as a measure of the physical size 
of the robot. Thus confining the robots to the e-robust 
visibility set guarantees free movement of the robot in 
the environment. Indeed, the notion of e-contraction is 
related to the classical work on motion planning in [15], 
see also [9]. □ 

We now define some graphs which will be useful in de- 
scribing the interactions between robots. 

Definition II.5 (Proximity graphs) A proximity graph is a 
graph whose nodes are a set of points V = {pi, . . . and 
whose edges are a function of V. Given V d Q, e > and 
r > 0, define: 

(i) The visibility graph ^/vis.Q at V is the graph with node 
set V and with edge set £vis..q{'P) defined by: {pi,pj) G 
^Vis^qiV) if and only f [pi,Pj] C Q. 

(ii) The e-robust visibility graph 5vis,Q, is the visibility graph 
at V for Qj. 

(iii) The r-disk graph Qr-a^k at V is the graph with node set 
V and with edge set fr-disk(7') defined by: {pi,Pj) G 
^r-diskiV) if and only if \\pi - pj\\ < r. 

(iv) The r-range visibility graph Qr-ma, q at V is the graph 
with node set V and with edge set £'r-disk('P) n^vis.gC'P)- 

(v) A Euclidean Minimum Spanning Tree t/EMST.g at V of 
a proximity graph Q is a minimum-length spanning tree 
of Q{V) whose edge (j)i,pj) has length \\pi — Pj\\. If 
Q{V) is not connected, then ^EMST,e(^) the union 
of Euclidean Minimum Spanning Trees of its connected 
components. 

In other words, two points p, q are neighbors in the r-range 
visibility graph, for instance, if and only if they are mutually 
visible and separated by a distance less than or equal to r. 
Example graphs are shown in Figure |4] General properties of 
proximity graphs are defined in [14], [7]. For simplicity, when 
Q is the complete graph, we denote the Euclidean Minimum 
Spanning Tree of Q by Gemst- 

We say that two proximity graphs Qi and Q2 have the 
same connected components if, for all sets of points V, the 
graphs Gi{V) and G2{'P) have the same number of connected 
components consisting of the same vertices. 

Definition II.6 (Neighbors set) Given a set of points V ~ 
{pi,...,p„} and a proximity graph Q, we let Mi{Q^V) ~ 
Mi.g{V) denote the set of neighbors including itself of pi. In 
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Other words, if {pi^ . 
V, then M{g,V) 



■ • ■ iPim} '^''^ ^''^ neighbors of pi in Q at 

M,e(P) = {p^„...,p^,JU{p^}. 





Fig. 4. The figure on tlie left stiows tiie visibility graph (whose edges are 
the solid lines as well as the dashed lines) and the e-robust visibility graph 
(whose edges are the solid lines alone) of a set of points in a nonconvex 
polygon. The figure on the right shows the r-range e-robust visibility graph. 
The disk in the figure shows the sensing range for one of the agents. 



Definition II.7 (Relative convex liull) Take an allowable 
environment Q. 

(i) X (1 Q is relatively convex ;/ the shortest path inside 
Q connecting any two points of X is contained in X. 

(ii) The relative convex hull rco{X, Q) of X C Q is the 
smallest relatively convex subset of Q that contains X. 

(iii) If X is a finite set of points, then a vertex o/rco(X, Q) 
is a point p €: X with the property that tco{X \ {p}, Q) 
is a strict subset of rco(X, Q). The set of vertices of 
rco(X, Q) is denoted by Ve(rco(X, Q)). 

The relative convex hull of an example set of points and its 
vertices are shown in Figure |5] In what follows we will need 




Fig. 5. Relative convex hull rco(X, Qe) of a set of points X (solid 
disks) inside a the e-contraction of an allowable set Q. The set of vertices 
Ve(rco(X, Qt)) is the set {v\ , . . . , fy}. 

the notion of perimeter of certain sets, and in particular, of the 
relative convex hull of a collection of points. 

Definition II.8 (Perimeter) Take an allowable environment 
Q and a closed subset X C Q. 

(i) If X has measurable boundary dX and is equal to the 
closure of its interior, then pcrimctcr(X) is the length 
ofdX. 

''That is, rco{X, Q) is the intersection of all relatively convex subsets of 
Q that contain X. 



(ii) If vco{X, Q) is not equal to the closure of its interior, 
then pcrimctcr(rco(X, Q)) is the length of the shortest 
measurable curve inside Q enclosing X. 

Remarks II.9 (i) If rco(X, Q) is equal to the closure of 
its interior, then its boundary is the shortest measurable 
curve inside Q enclosing X (i.e., the two definitions 
of perimeter are equivalent). On the other hand, if 
rco(X, Q) is a segment, then Definition III. 8f ii) says that 
the perimeter of rco(X, Q) is twice its length, 
(ii) The key property of Definition III. 81 is that, if X is a 
finite set of points in Q, then the perimeter of rco(X, Q) 
depends continuously on the points m X. □ 

III. Synchronous robots with visibility sensors 

AND THE rendezvous AND CONNECTIVITY 
MAINTENANCE PROBLEMS 

In this section we model a group of n robots with visibility 
sensors in a given allowable environment Q. We assume that 
e is a known positive constant sufficiently small so that Q^^ 
is allowable. For i € {!,..., n}, we model the ith robot as a 
point Pi e Q and we refer to Section |VT] for an extension to 
a disk model. We make the following modeling assumptions: 
Synchronized controlled motion model: Robot i moves at 
time t S Z>o for a unit period of time, according to 
the discrete-time control system 



p^[t + I] ^ Pi[t] + u^[t\. 



(1) 



We assume that there is a maximum step size s^ax > for 
any robot, that is, \\ui\\ < Smax- Note that the n identical 
robots are synchronized in the sense that the calculation 
of u[t] in equation ([T]) takes place at the same times t for 
all robots. 

Sensing model: Robot i senses (i) the presence and the 
position of any other robot that is visible and within 
distance r from pi, and (ii) the subset of dQ that is 
visible and within distance (r + e) from pi. This in turn 
implies that the robot can sense the subset of dQ^ that is 
visible and within distance r from pi. It is convenient to 
define the sensing region from position pi to be S{pi) = 
V{pi, e)r\B{pi, r). The range r is the same for all robots. 
Note that, by definition, two robots with visibility sensors 
detect each other's presence and relative position if and only 
if they are neighbors in the robust visibility graph ^vi,s,Q5 ■ 

Remark III.l (No common reference frame) The model 
presented above assumes the ability of robots to sense 
absolute positions of other robots; this assumption is only 
made to keep the presentation as simple as possible. In this 
and subsequent remarks, we treat the more realistic setting 
in which the n robots have n distinct reference frames 
I]i,...,S„. We let Eo denote a fixed reference frame. 
Notation-wise, a point q, a vector w, and a set of points S 
expressed with respect to frame are denoted by g*, 
and 5*, respectively. For example, this means that is the 
environment Q as expressed in frame E^. We assume that the 
origin of E; is pi and that the orientation of E; with respect 
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to Eq is i?" G SO(2). Therefore, changes of reference frames 
are described by the equations: = i?,"(7' +P^, w° = R^w\ 
and S"" = i?,"S'* If we let Vqi (p^ , e) denote the visibiUty 
set expressed in Ej, for G {0, 1, . . . , ??}, then one can define 

S{plQ^)^VQ,{ple)r\B{plr), 

and verify that Q") = Q')+p°. Note that = 

0. 

Finally, we can describe our motion and sensing model 
under the no common reference frame assumption. Robot i 
moves according to 

pO[t+l]=p^[t] + i?,"[i]u,[i], (2) 

and it senses the robot positions and the subset of {dQY 
that are within the sensing region S{pl, Q*). □ 

We now state the two control design problems addressed in 
this paper for groups of robots with visibility sensors. 

Problem III.2 (Rendezvous) The rendezvous problem is to 
steer each agent to a common location inside the environ- 
ment Qj. This objective is to be achieved (1) with the limited 
information flow described in the model above, and (2) under 
the reasonable assumption that the initial position of the robots 
P[Q] = {pi[0], . . . ,p„[0]} gives rise to a connected robust 
visibiUty graph 5vis,Q, at V[0]. □ 

As one might imagine, the approach to solving the ren- 
dezvous problem involves two main ideas: first, the underlying 
proximity graph should not lose connectivity during the evo- 
lution of the group; second, while preserving the connectivity 
of the graph, the agents must move closer to each other. This 
discussion motivates a second complementary objective. 

Problem III.3 (Connectivity maintenance) The connectiv- 
ity maintenance problem is to design (state dependent) control 
constraints sets with the following property: if each agent's 
control takes values in the control constraint set, then the 
agents move in such a way that the number of connected 
components of tJvis.g, (evaluated at the agents' states) does 
not increase with time. □ 

IV. The connectivity maintenance problem 

In this section, we maintain the connectivity of the group of 
agents with visibility sensors by designing control constraint 
sets that guarantee that every edge of Qr-via, (i-S-, every pair 
of mutually range-limited visible robots) is preserved. We have 
three objectives in doing so. First, the sets need to depend 
continuously on the position of the robots. Second, the sets 
need to be computed in a distributed way based only on the 
available sensory information. Third, the control constraint sets 
should be as "large" as possible so as to minimally constrain 
the motion of the robots. Because it appears difficult to 
formalize the notion of "largest continuous constraint set that 
can be computed in a distributed fashion," we instead propose 
a geometric strategy to compute appropriate constraint sets 
and we show in the next section that our proposed geometric 
strategy is sufficiently efficient for the rendezvous problem. 



A. Preserving mutual visibility: The Constraint Set 
Generator Algorithm 

Consider a pair of robots in an environment Q that are e- 
robustly visible to each other and separated by a distance not 
larger than r. To preserve this range-limited mutual visibility 
property, we restrict the motion of the robots to an appropriate 
subset of the environment. This idea is inspired by [5] and we 
begin by stating the result therein. Let the sensing region of 
robot i located at pi be S{pi) — B{pi, r), for some r > 0. If 
at any time instant t, \\pi[t] — Pj[t]\\ < r then to ensure that 
at the next time instant t + 1, \\pi[t + 1] — Pj[t + 1]|| < r, it 
suffices to impose the following constraints on the motion of 
robots i and j: 

or, equivalently, 

In summary, B{ ^^~^^ , |) is the control constraint set for robot 
i and j. This constraint is illustrated in Figure |6] (left). 




Fig. 6. In the figure on the left, starting from pi and pj, the robots 
are restricted to move inside the disk centered at ^^iih radius 

I . In the figure on the right, the robots are constrained to move inside 
the shaded region which is a convex subset of intersected with 
the disk centered at ^i-tZi ^^Jth radius §. 

Let us now consider the case when a robot i is located at pi 
in a nonconvex environment Q with sensing region S{pi) = 
V{pi, e)r\B{pi, r). If at any time instant t, we have that 
Pj[t]\\ < r and G Qe, then to ensure that + 

1] - pj[t + < r and [p^[t + l],pj[t + 1]] G Qe, it suffices 
to require that: 

P,[t + 1], pj[t + i]eC, 

where C is any convex subset of n ^^(^PiMliLLW ^ i^; see 
Figure |6] (right). Equivalently, 

u,[t] eC^p.,[t], uj[t] eC-pj[t], 

where C — pi[t] and C — pi[t] are the sets {p — Pi[t] | p G C} 
and {p — pj[t] I p G C}, respectively. Note that both robots 
i and j must independently compute the same set C. Given 
the positions pi , pj in an environment Q, Table U describes 
the Constraint Set Generator Algorithm, a geo- 
metric strategy for each robot to compute a constraint set 
C = CQ{pi,pj) that changes continuously with pi and pj. 
Figure |7] illustrates a step-by-step execution of the algorithm. 



TABLE I 

Constraint Set Generator Algorithm 

Goal: Generate convex sets to act as constraints to preserve mutual 
visibility 

Given: (Pi,Pj) G such that [pi,Pj] C and pj e B{pi,r) 
Robot i £ {1, . . . ,n} executes the following computations: 

1: C.emp := V(pi,e)nB(^^,§) 

2: while 9Ctcmp contains a concavity do 

3: ti : = a strictly concave point of cJCtcmp closest to the segment [pi,Pj] 

4: Ciemp := Ctcmp H Hq^ (v) 

5: end while 

6; return: CQ(pi,pj) := Ctcmp 




Fig. 7. From left to right and top to bottom, a sample incomplete run of the 
Constraint Set Generator Algorithm (cf. Table H). The top left 
figure shows Ctcmp := V(pi,e) n , ^). In all the other figures, the 

lightly and darkly shaded regions together represent Ctcmp. The darkly shaded 
region represents Ctcmp H Hq(v), where v is as described in step 3 : . The 
final outcome of the algorithm, CQ{pi, pj), is shown in Figure |6] (right). 

In step 3 : of the algorithm, note that there can be multiple 
distinct points belonging to distinct concavities satisfying the 
required property. In that case, v can be chosen to be any one 
of them. The following lemma justifies this observation. 

Lemma IV.l Throughout the execution of the Constraint 
Set Generator Algorithm in Tabled let Vi, V2 be two 
strictly concave points on 9Ctemp that are closest to [pi,pj\. 
Then Vi £ Ctemp H Hq^ (^2) and vice versa. 

Next, we characterize the main properties of the 
Constraint Set Generator Algorithm and 
the corresponding convex sets. Notice that the constraint set 
is defined at any point in the following set: 

J = {{pi,Pj) e Ql I [pi,Pj] e Qe, \\pt -P]\\ < r}. 



Proposition rv.2 (Properties of the Constraint Set 
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Fig. 8. The green convex set in the center represents Cp^^Q{N'i^g^^^^ ^ )■ 
The black disks represent the position of the robots. The straight line segments 
between pairs of robots represent edges of 5r-vis,Q ■ Here, pi is the black disk 
contained in the constraint set. 

Generator Algorithm) Given an allowable environ- 
ment Q with K strict concavities, e > and {pi,Pj) G J, 
the following statements hold: 

(i) The Constraint Set Generator Algorithm 
terminates in at most k steps; 

(ii) Cq {pi , Pj ) is nonempty, compact and convex; 

(iii) CQ{pi,pj) =CQ{pj,pi); and 

(iv) The set-valued map Cq is closeiQ at every point of J. 

Remark IV.3 (No common reference frame: continued) 

Consider a group of robots with visibility sensors and no 
common reference frame. With the notation and assumptions 
described in Remark IIII.U on can verify that the constraint 
sets transform under changes of coordinate frames according 
to; 

CQo(pO,pO) = i^^CQ.(p^pj)+rf. (3) 

We omit the proof in the interest of brevity. □ 

For each pair of mutually visible robots, the execution of 
the Constraint Set Generator Algorithm outputs 
a control constraint set such that, if the robots' motions are 
constrained to it, then the robots remain mutually visible. 
Clearly, given a connected graph at time t, if every robot 
remains connected with all its neighbors at time t + 1 (i.e., 
each pair of mutually visible robots remain mutually visible), 
then the connectivity of the graph is preserved. This can be 
accomplished as follows. For robot i G {1, . . . , n} at pi G Q^, 
define the control constraint set 

C^P.,g(M.e...,Qj = n Cq(p„p,). (4) 

Now, if Ui e C'p,,Q(A/'i,e.-vis. Q, ) -Pi^ for alH e {1, . . . , n}, 
then all neighboring relationships in Qr-wis, ^re preserved 
at the next time instant. Using inputs that satisfy these con- 
straints, the number of edges in Gr-vis, is guaranteed to be 
nondecreasing. 

^Let fl map points in X to all possible subsets of Y. Then the set- valued 
map, f2, is open at a point x G X ii for any sequence {x^} in X, — > x 
and y S Q{x) implies the existence of a number m and a sequence {yk} 
in Y such that yf^ g ^(^k) for k > m and — > y. The map f! is closed 
at a point x £ X ii for any sequence {x't} in X, Xf^ — ► x, yk ^ y and 
Vk G ^{^k) imply that y £ f2(x). f7 is continuous at any point a; S X if it 
is both open and closed at x 
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B. The locally-cliqueless visibility graph 

In this section, we propose the construction of constraint 
sets that are, in general, larger than Cp-^Q{N'i.g^_^-^ q^). To 
do this, we define the notion of locally-cliqueless graph. The 
locally-cliqueless graph of a proximity graph Q is a subgraph 
of Q, and therefore has generally fewer edges, but it has the 
same number of connected component as Q. This fundamental 
property will be very useful in the design of less conservative 
constraint sets. 

Before proceeding with the definition of locally-cliqueless 
graph, let us recall that (i) a clique of a graph is a complete 
subgraph of it, and (ii) a maximal clique of an edge is a clique 
of the graph that contains the edge and is not a strict subgraph 
of any other clique of the graph that also contains the edge. 

Definition IV.4 (Locally-cliqueless graph of a proximity 
graph) Given a point set V and an allowable environment Q, 
the locally-cliqueless graph Q\c^g at V of a proximity graph Q is 
the proximity graph with node set V and with edge set f ic,c; (7-") 
defined by: {pi,Pj) G Sgj^.gCP) if and only if {pi,Pj) G ^giV) 
and (pi , pj ) belongs to a set Sg^f^si i'P') fi"" '^^y maximal clique 
V' of the edge {pi , pj ) in Q. 



In combinatorial optimization, it is a well-known that find- 
ing the maximal clique of a graph is an NP complete problem. 
However, efficient polynomial time heuristics are detailed 
in [16]. 

For simplicity, we will refer to the locally-cliqueless graph 
of the proximity graphs t/vis,Q, 0vis,Q. or t/r-vis,Q, as locally- 
cliqueless visibility graphs. Figure |9] shows an example of a 
locally-cliqueless visibility graph. 




Fig. 9. Visibility grapli (left) and locally-cliqueless visibility graph (right). 



Theorem IV.5 (Properties of a locally-cliqueless graph of 
a proximity graph) Let Q be a proximity graph. Then, the 
following statements hold: 

(i) Gemsw ^ S\c,g C G; 

(ii) Gic.g ond G have the same connected components. 

In general, the inclusions in Theorem |IV.5[ i) are strict. Fig- 
ure [To] shows an example where ^/emst.s.^.q £ Gic,g,„,Q Q 

Gms.Q- 

The next result follows directly from Theorem IIV.5I 

Corollary IV.6 (Properties of locally-cUqueless visibility 
graphs) Let Q and Q^, e > 0, be allowable environments. 
Let G be either one of the graphs Gyis,Q, Gns,Qe or Sr-vis.Q^- 
Then, the following statements hold: 




Fig. 10. From left to right, visibility graph, locally-cliqueless graph and 
Euclidean Minimum Spanning Tree of the visibility graph. 



(i) Semst.s C Q\^ g C Q; 

(ii) Gxc.g cind G have the same connected components. 

Let us now proceed to define new constraint sets that are in 
general larger than the ones defined in (|4|. For simplicity, let 
G = ^r-vis.Q^, and consider its locally-cliqueless graph G\c,g- 
For robot i G {1, . . . , n} at position pi, define the constraint 
set 

Cp.,Q(M,e,c,s)= n Cq(p„p,). (5) 

Since Qx^^ g is a subgraph of Q (cf. Corollary IIV.6f i)), we have 



and therefore 



In general, since Qx^^^g is a strict subgraph of Q, the set 
Cpi,g(-^i,ek.s) is strictly larger that Cp,,g(A/'i,g^_„, ^J. 

Now, if Ui g Cp.,Q(A/i,g„ g)— Pi for alH S {1, . . . , n}, then 
all neighboring relationships in the graph t/ic, g are preserved 
at the next time instant. As a consequence, it follows from 
Corollary IIV.6t ii) that the connected components of Gr-ma, 
are also preserved at the next time instant. Thus, we have 
found constraint sets (|5]) for the input that are larger than 
the constraint sets (|4]i, and are yet sufficient to preserve the 
connectivity of the overall group. 

Remark IV.7 (Distributed computation of locally- 
cliqueless visibility graphs) According to the model 
specified in Section each robot can detect all other robots 
in its sensing region S{pi) = V(pi,e) fl B{pi,r), i.e., its 
neighbors in the graph Gr~vis, q, ■ Given the construction of the 
constraint sets in this section, it is important to guarantee that 
the set of neighbors of robot i in the locally-cliqueless graph 
Qxi^^g can be computed locally by robot t. From the definition 
of the locally-cliqueless graph, this is indeed possible if a 
robot i can detect whether another robot j in its sensing 
region S{pi) belongs to a clique of the graph Gr-ms, q,- This is 
equivalent to being able to check if two robots putPi G S{pi) 
satisfy the condition that pk G S{pi) and vice versa. Note that 
Pk G S{pi) is equivalent to \\pk - pi\\ < r and [pk,Pi] G Q^. 
Given that Pk - Pi = {pk - Pi) - (pi -Pi), the vector pk - pi 
(and hence \\pk — pi\\) can be computed based on local 
sensing alone. Now, checking if [pk,Pi] C is possible 
only if does not contain any hole; see Figure [TT] In such 
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TABLE II 

Perimeter Minimizing Algorithm 



Fig. 11. The dashed circle is centered at pi and is of radius r. The thick 
curves represent the boundary of Q^', the one on the left represents the outer 
boundary whereas the one on the right represents a hole in the environment. 



a case, it suffices to check if the entire line segment [pk,pi] 
is visible from pi or not. 

Along these same lines it is possible to state that the locally- 
cliqueless visibility graph can be computed also under the "no 
common reference frame" model described in Remarks IIII.ll 
andHVj] □ 

V. The rendezvous problem: Algorithm design and 

ANALYSIS RESULTS 

In this section, we solve the rendezvous problem through a 
novel Perimeter Minimizing Algorithm. The algo- 
rithm is inspired by the one introduced in [5] but is unique in 
many different ways. The rendezvous algorithm uses different 
graphs to maintain connectivity and to move closer to other 
robots. Instead of moving towards the circumcenter of the 
neighboring robots, the robots move towards the center of a 
suitably defined motion constraint set. 

The section is organized as follows. We present the algo- 
rithm in Subsection IV-AI followed by its main convergence 
properties in Subsection IV-BI 

A. The Perimeter Minimizing Algorithm 

We begin with an informal description of the Perimeter 
Minimizing Algorithm over graphs Qsens and ^/constr- 
The sensing graph Q^ens is Sr-vis.g, while the constraint graph 
Sconstr is either g^tns or Gicg,,^,'- 

Every robot i performs the following tasks: (i) it 
acquires the positions of other robots that are its 
neighbors according to Q^ens', (ii) it computes a 
point that is "closer" to the robots it senses, and 
(iii) it moves toward this point while maintaining 
connectivity with its neighbors according to ^?consti - 
The algorithm is formally described in Table [III Figure [1] in 
the Introduction illustrates an example execution. 

Remarks V.l (i) According to the algorithm proposed 
in [5] the robots move towards the circumcenter of their 
neighbors position. In the Perimeter Minimizing 
Algorithm the robots move towards the circumcenter 
of their constraint set. 
(ii) We prove later in Lemma HOI in Appendix HIl that Xi is 
convex and that, in turn, CC{Xi) is well-defined. Be- 
cause GC{Xi) S Xi, we know that p* € Xi. Therefore, 
Ui € Xi - pi C Cp^,Q{Aft,g,„„,^,) - Pi and, in turn. 
Pi at the next time instant belongs to C'p..Q(A/i.g„„,J. 



Assumes: (i) 

(ii) 



Smax > is the maximal step size 
Q,Qc are allowable 
(iii) Gxm is gr-vis. Qr '■ gconstr is either gsens or Gic,g„,„: 



Each robot i G {1, . . . , n} executes the following steps at each time instant: 
1: acquire {pi-^ , ■ ■ ■ , Pim } ~ positions of robots within pi sensing region 



compute Ni^g^^, and M.Sconar 

compute Xi := Cp;^Q(A/'i,g„„„J n rco(A/'i,6;,„, , V(pi, €)) 
compute p* := CC{Xi) 

_ min(5max, lb* -PiW) , 

M - Pill 



2 
3 
4: 

5: return: Ui 



-(P*i - Pi) 



From our discussion in Section IIVI this implies that the 
graph ^constr remains connected (or, more generally, that 
the number of connected components of ^constr does not 
decrease). Therefore, by Corollary IIV.6I the number of 
connected components of C/sens also does not decrease. 

(iii) If the initial positions of the robots are in (Je, then the 
robots will remain forever in because p* & Xi <Z Q^. 

(iv) All information required to execute the steps in the 
algorithm is available to a robot through the sensing 
model described in Section [Till The constraint on the 
input size, < Smax, is enforced in step 5:. □ 

Finally, we conclude this section by completing our treat- 
ment of robots without a common reference frame. 

Remark V.2 (No common reference frame: continued) 

Consider a group of robots with visibility sensors and no 
common reference frame as discussed in Remarks IIII.ll 
and IIV.3I Because the relative convex hull and the 
circumcenter of a set transform under changes of coordinate 
frames in the same way as the constraint set does in 
equation (O, one can verify that 

where Ui{p1, . . . is computed with environment Q'^ and 
Ui{p\, . . . ,pU is computed with environment Q^. This equal- 
ity implies that the robot motion with control Ui{j)\, . . . ,p°) 
in equation ([T]i is identical to the robot motion with control 
Ui{p\, . . . in equation □ 

B. Main convergence result 

To state the main results on the correctness of the 
Perimeter Minimizing Algorithm, we require some 
preliminary notation. First, note that given the positions of 
the robots {pi, . . . ,p„} at any time instant t, the algorithm 
computes the positions at time instant i + 1 . We can therefore 
think of the Perimeter Minimizing Algorithm as 
the map Tg ^.^ ^g^^,^^,^ : (J" Q^. Second, in what follows we 
will work with tuple of points P = {pi, . . . ,pn) S Q"- We let 
Q{P) denote the proximity graph 0(1^) and rco(P, Q) denote 
the relative convex hull of the set V inside Q, where V is the 
point set given by {pi \ i G {l,...,n}}. Third, we introduce 
a Lyapunov function that encodes the rendezvous objective. 
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Given an allowable environment Q, we recall the notions of 
relative convex hull and of perimeter from Section |II] and we 
define Vperim.g : Q" ^ K>o by 

Vperim,Q(-P) = perimeter(rco(P, Q) ) . 



Lemma V.3 (Properties of Lyapunov function) The func- 
tion Vjjerim.Q has the following properties: 

(i) V^erim.Q is continuous and invariant under permutations 
of its arguments; 

(ii) Vperim,g(P') = /or P = (pi , . . . , p„) if and only if 
Pi = pj for all i, j E {1, . . . , n}. 

The technical result on continuity plays an important role 
in the convergence analysis of the algorithm. Fact (ii) implies 
that achieving the rendezvous objective is equivalent to making 
V^erim.Q, equal to zero. We are now ready to state the main 
result of the paper. 

Theorem V.4 (Rendezvous is acliieved via tlie Perimeter 
Minimizing Algorithm) Let Q and Qf be allowable en- 
vironments. Let pi , . . . , p„ be a group of robots with visibility 
sensors in Q^. Any trajectory C Qe generated by 

P[t + 1] = Tescns.ec™sir(-PM) ^'^^^ thc followlng propcrtlcs: 

(i) ;/ the locations of two robots belong to the same con- 
nected component of Q^ens cit P[tQ\ for some to, then 
they remain in the same connected component of Qaena 
at P\t] for all t > to," 

(ii) ^perim,Q.(-P[i + 1]) < '^'perim.Q. (^^M )/ and 

(iii) the trajectory {P[t]}tei.>o converges to a point P* £ Qe 
such that either p* = p* or p* ^ S{Pj) for all 
{!,..., n}. 

As a direct consequence of the theorem, note that if the graph 
Gsens is connected at any time during the evolution of the 
system, then all the robots converge to the same location in Q^. 

VI. Practical implementation issues 

In Section |V] we designed a provably correct rendezvous 
algorithm for an ideal model of point robots with perfect 
sensing and actuation capabilities. Also, we assumed that it 
was possible for the robots to operate synchronously. However, 
such an ideal model is not realistic in practical situations. In 
this section, we investigate, via extensive computer simula- 
tions, the effects of deviations from this ideal scenario. 

A. Nominal experimental set-up 

The computer simulation was written in C++ using the 
Computational Geometry Algorithmic Library (CGAL) 
( |http : / /www . cgal . org| . However, it was found 
that Boolean operations on polygons using the utilities 
present in CGAL were not adequate in terms of 
speed for the purpose of running extensive simulations. 
Hence, Boolean operations on polygons were performed 
using the General Polygon Clipping Library (GPC) 
( |http : / /www . cs ■ man . ac . uk/~toby7alan/ software 



For the purpose of simulations, the environment considered 
is a typical floor plan; see Figure [T] Experiments were per- 
formed with 20 robots starting from 10 randomly generated 
initial conditions from a uniform distribution. For each initial 
condition, experiments were repeated 20 times. The environ- 
ment size is roughly 80 x 70, the step size of a robot is taken 
as Smax = 0.5 and the sensing radius r = 30. For simplicity, 
Sconstr is taken to be the same as Ssens- To utilize the e-robust 
visibility notion in providing robustness to asynchronism and 
sensing and control errors, at each time instant, e is set to 
be 0.97 times the e at the previous time instant. Initially, 
e is set equal to 3. In case a robot approaches a reflex 
vertex of the environment closely, it reduces its speed. This is 
done to reduce the risk of collision due to errors in sensing 
the exact location of the reflex vertex. We now describe 
the various assumptions we make on the model to simulate 
the actual implementation on physical robots followed by 
the respective simulation results. The algorithm performance 
is then evaluated based on the following three performance 
measures: (i) the average number of steps taken by the robots 
to achieve the rendezvous objective; (ii) the fraction of the 
edges of C/sens that are preserved at the end of the simulation; 
and (iii) the number of connected components of CJsens at the 
end of the simulation compared with the number of connected 
components at the initial time. 

B. Robustness against asynchronism, sensing and control 
noise, and finite-size disk robot models 

(Al) Asynchronism: The robots operate asynchronously, i.e., 
do not share a common processor clock. All the robots 
start operating at the same time. Each robot's clock 
speed is a random number uniformly distributed on 
the interval [0.9,1]. At integral multiples of its clock 
speed, a robot wakes up, senses the positions of other 
robots within its sensing region and takes a step accord- 
ing to the Perimeter Minimizing Algorithm. 
Note that at the time when any given robot wakes up, 
there may be other robots that are moving. No sensing 
and control errors were introduced in the model. It is 
observed that under this asynchronous implementation, 
the performance of the algorithm is very similar to 
the synchronous implementation; see Figure [12] In the 
subsequent implementations, we assume the robots to 
operate synchronously. 

(A2) Distance error in sensing and control: The visibility 
sensors measure the relative distance of another object 
according to the following multiplicative noise model. 
If dact is the actual distance to the object, then the 
measured distance is given by (1 + edist)rfact, where 
edist is a random variable uniformly distributed in the 
interval [—0.1,0.1]. The objects to which distances are 
measured are other robots in the sensing range and the 
boundary of the environment. For simplicity, instead of 
measuring a sequence of points along the boundary, as a 
real range sensor does, we assume that only the vertices 
of the environment are measured and the sensed region 

Tgpc ijhtacclt| listructed from that information. The sensor error. 
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Fig. 12. Computer simulation results with asynchronism. The top figure rep- 
resents the average number of steps taken per robot for convergence(crosses). 
Also shown is the number of steps taken by each robot in synchronous 
implementation (red circle). The center figure shows the fraction of the 
edges of the initial sensing graph that are preserved till the end. The bottom 
figure shows the number of connected components of the sensing graph 
initially (small black discs) and at the end of asynchronous (blue crosses) 
and synchronous (red circle) implementations. The blue crosses in the figure 
denote the mean of the observed quantities and the vertical bars denote 
standard deviation. 



therefore, occurs in the measurement of other robots and 
environment vertices. The actuators moving the robots 
are also subject to a muhipHcative noise distance model 
with the same error parameter The results are shown 
in Figure [13] It is observed that only in 1 out of 



identical to the synchronized implementation with no 
noise. In the subsequent simulations, we assume the 
robots to operate synchronously with no distance error 
in sensing and control. 
(A3) Direction error in sensing and control: The visibility 
sensors measure the relative angular location of another 
object according to the following additive noise model. 
If 0act is the actual angular location of any object in the 
local reference frame of a robot, the measured angular 
location is given by ^act + ee, where eg is a random 
variable uniformly distributed in the interval [—5, 5]. As 
before, the actuators moving the robots are also subject 
to an additive noise directional model with the same 
error parameter The results are shown in Figure [l4l 
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Fig. 13. Computer simulation results with distance error and no directi onal 
error in sensing and control. The meaning of the quantities is as in Figure 1121 



the 10 initial conditions does the number of connected 
components of tJsens increase as compared to the number 
of connected components of t/jens initially. In all cases, 
the fraction of edges preserved is almost equal to 1. 
Also, in 7 out of 10 cases, the performance is almost 



Fig. 14. Computer simulation results with direction en'or and no distance 
en'or in sensing and control. The meaning of the quantities is as in Figure [12] 

It is observed that the algorithm no longer performs 
similar to the synchronized implementation with no 
noise. Thus, it is more sensitive to directional errors 
than distance errors in sensing and control. However, 
almost all of the edges of Gsens are preserved for all 
the initial conditions. Also, in n 9 out of the 10 initial 
conditions the number of connected components of (?sens 
corresponding decreases during the course of evolution 
of the group. In the following simulations, no asynchro- 
nism or sensing and control errors are assumed. 
(A4) Disk robot model: The robots are assumed to be disks 
of radius 0.2, but do not obstruct the views of others. 
During any step, a robot moves a distance of at most 
Smax = 0.5. The next position of the center of the 
robot, therefore, lies in a motion disc of radius 0.7 
centered at the center of the robot; see the red and 
green discs in Figure [15] A colliding neighbor of a 
robot i is any neighbor j according to f/sens such that 
the motion discs of i and j intersect and that the motion 
disc of j intersects the physical disc of i on the path 
between i and its next point. If a robot has no colliding 
neighbors, then its motion is according to Perimeter 
Minimizing Algorithm. If on the other hand a 
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robot has exactly one colliding neighbor, then it tries 
to swerve around it while reducing its speed. Finally, if 
the robot has more than two colliding neighbors then it 
stays at the current location. To ensure free movement of 
the robots inside the environment, e is not allowed to fall 
below 0.2, i.e., the radius of a robot disc. Simulations 
were performed without any asynchronism or sensing 
and control errors for the 10 initial conditions of the 
robots as in the previous experiments. The terminating 
condition for the simulations is that robots belonging to 
each connected component of C/sens form a "cohesive" 
grouf@; see Figure [15] For all initial conditions, the 
number of cohesive groups in the final configurations 
is equal to the number of connected components of 
Gsens when the simulations are performed assuming point 
robots. Thus, the algorithm yields the same performance 
if a disk robot model is assumed instead of a point robot 
model. 




"Til 



Fig. 15. Computer simulation results with no asynchronism and no sensing 
and control eiTors. The black and red discs denote the robots and their motion 
discs respectively. Th e initia l p ositio n of the robots correspond to initial 
condition 2 in Figures [Til 1 13| and l 14l and are shown by the small black discs 
scattered over the environment with the green discs denoting their motion 
discs. The robots converge to positions corresponding to a single cohesive 
group part of the same component of Gsens- 

Thus, we see that the Perimeter Minimizing 
Algorithm is robust to various deviations from the 
ideal scenario. The magnitudes of edist and ee are in line with 
the state-of-the-art. Finally, in the next section we analyze 
the computation complexity of the algorithm. 

C. Computation complexity with finite resolution sensing 

In addition to the issues that might arise in a practical imple- 
mentation of the Perimeter Minimizing Algorithm, 
another important consideration is the time taken for a robot 
to complete each step of the algorithm. This is dependent 

*For each connected component of Sscns. the graph having nodes as 
the robot locations and with an edge between two nodes whenever the 
corresponding motion discs of the robots intersect is connected 



on the computational complexity of the algorithm, that we 
characterize in the following. 

A real visibility sensor, e.g., a range scanner, will sense the 
position of other robots and the boundary of the environment 
with some finite resolution; in particular, the boundary of 
the sensing region will be described by a set of points. It is 
reasonable to assume that the cardinality of this set of points 
is bounded, say by M, for all robots, irrespective of the shape 
of the environment and the location of the robot in it. For 
example, if a laser range sensor is used to measure the distance 
to the boundary and a measurement is taken at intervals of one 
degree, then AI is equal to 360. 

Proposition VI.l (Computational complexity) Let Q be 

any allowable environment. Let M be the resolution of the 
visibility sensor located at any robot in Q. Then the following 
statements are true: 

(i) The computation complexity of the Constraxnt. Set 
Generator Algorithm is 0{kM); 

(ii) The computation complexity of the Perimeter 
Minimizing Algorithm is t{M) + 0{M^ logM); 

(iii) If t/constr ^ Gsens, then the computation complexity 
of the Perimeter Minimizing Algorithm is 
0(Af2 logAf), 

where t{M) is time taken for the computation o/TVi^g;,.^^,, given 
the set Mi^g^^^^ assuming lA/i^g^^^ | < M, and k is the number 
of strict concavities of Q. 

As discussed in Section IIV-BI if ^constr = Q\c,g^„,, then the 
computation of ^constr from Gsens can be performed using 
efficient polynomial time heuristics. The time t{M) above 
depends on the specific heuristic used. Thus, we see that the 
running time of each step of the Perimeter Minimizing 
Algorithm is polynomial in the number of data points 
obtained by the visibility sensor. 

This analysis helps us assess whether it is feasible to im- 
plement this algorithm on an actual robot without demanding 
an unreasonable computational power. 

VII. Conclusions 

In this paper, we present a provably correct discrete-time 
synchronous Perimeter Minimizing Algorithm al- 
gorithm for rendezvous of robots equipped with visibility 
sensors in a nonconvex environment. The algorithm builds 
on a novel solution to the connectivity maintenance problem 
also proposed in this paper. The performance of the algorithm 
under asynchronous operation of the robots, presence of noise 
in sensing and control, and nontrivial robot dimension is inves- 
tigated and found to be quite satisfactory. The computational 
complexity of the algorithm under the assumption of finite 
sensing resolution is also investigated. 
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Appendix I 
Proofs of results in SectionITvI 

Proof of Lemma \IV.1\ Let d = dist{vi, [pi,pj\) = 
dist(w2, [pi,P]]) and let L = {p eM? \ dist(p, [pi.Pj]) < d}. 
Then vi , V2 S dL; see Figure [16] We now prove our result by 
contradiction. Let vi ^ Hq^{v2). Then vi E dL \ Hq^{v2) as 
shown in Figure [161 Since vi is visible from pi, the boundary 
dQe must intersect dL at a point v in between V2 and vi, 
as illustrated in Figure [16] But this means that there exists 
a point V* belonging to the concavity containing V2 that 
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Fig. 16. The shaded region represents L. The solid curve passing through V2 
represents a portion of the boundary of Qe . The dashed line is the boundary of 
Hq^ {v2) which is along the tangent to dQ^ at V2. The interior of Hq^ (112) 
is in the direction of the arrow. 

is Strictly closer to the segment than V2', this is a 

contradiction. ■ 
Proof of Proposition \IV.2\ We first prove statement (i). Note 
that initially Ctemp = V(pi, e) n B{^^, §). Therefore (9C,emp 
has at most k concavities. This is because the concavities of 
9Ctemp are induced by the concavities of dQe, and the number 
of concavities of dQe is the same as the number of concavities 
of dQ. Now, by construction, the number of concavities in 
(?(Ctemp n Hq^ (v)) is strictly less than the number of concav- 
ities of dCtemp- It follows then that the Constraint Set 
Generator Algorithm terminates in at most k steps. 
This completes the proof of statement (i). 

Now note that [pi,Pj] C Cq{pi,Pj). Hence Cq{pi,Pj) is 
always non-empty. Notice that CQ(pi,pj) can be written in 
the following way 

CQ(p„p,) = V(p„6)nB(5i±^,^) 

n Hq, ...nHg, {■V^,J,n,^), (6) 

where dist(uij,i, [pj,pj]) < ... < dist(ujj_„,^ , 
The Wi,j.fc's are computed according to step 3: in the 
Constraint Set Generator Algorithm. Therefore, 
CqiPiiPj) is the result of the intersection of a finite num- 
ber of closed sets, and therefore is closed as well. Also 
CqiPiiPj) C and is thus bounded. Therefore, it is compact. 
Finally, the Constraint Set Generator Algorithm 
terminates when Ctemp has no concavities, or in other words, 
when Ctemp is convex. This proves statement (ii). 
To prove statement (iii), let us write CQ{pj,pi) as 

CQ(p„pO = V(p.,e)nB(5i±^,^) 

in the same way as We first claim that riij = nji and 
{vi^j^i, . . . ,Vi^j^ni,} = {vj,i,i, ■ ■ ■ ,Vj.i^nji}- Before proving 
this claim, let us assume it is true and see where it leads. 
Because of this assumption, the expressions for CQ{pi,Pj) 
and CQ{pj,Pi) would differ only due to the terms V{pi,e) 
and V{pj,e). Now, let q S CQ{pi,Pj). Because Cqip.i^Pj) is 
a convex subset of containing pj and q, the line segment 
[pj,q] C CQ{pi,Pj) C Q^. Hence, q £ V{pj,e). This in turn 
implies that q S CQ{pj,pi). Therefore, we have CQ{pi,pj) C 
CQ{pj,pi). The opposite inclusion can be shown to be true in 
a similar fashion. Thus, we have CQ{pi,pj) = CQ{pj,pi). 

We now prove that = nji and . . . ,Wi.j,Tii } = 

{vj^i^i, . . . jVj^i^n-i}- Note that is the strictly concave 

point nearest to [pi,Pj] belonging to V{pi,e) D ^ r). 

In general, there can be more than one nearest strictly concave 
point, say . . . Let d = dist(wij4, [PiyPj]) 
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L = {p' e CtempI dist(p', [pj,Pj]) < d], where Ctemp = 
V{pi,e) n , j). Note that L is convex since there 

cannot be any strictly concave points of i9Ctemp in the interior 
of L. Since pj, i^ij i G L then [vij,i^pj\ C L C Q^. Therefore 
w^.j-i e V(pj,e) and hence e V(pj,e) nB(^^4^,5). 

Therefore, u.j- 1, . . . , t.^^- £ V(pj, e) n §). Sim- 

ilarly, if Vj^i^i, . . . ,Vj^i^sji are the strictly concave points 
nearest to [pj,Pj] belonging to V{pj,e) fl ^ r)^ 

then . . . , G V(pi, e) n §). Therefore 

dist(u,;j_i, [pi,pj]) = dist(wj,i,i, and = Sj,. 

This implies that Wj,i,fc £ {wi^.i, . . . , j .,. . } for all fc G 
{1, . . . , Sji}. Proceeding recursively, we get that ri.y = rij^ and 
. . . , = ■ • ■ , ^'j.^niJ- This completes 

the proof of statement (iii). 

We now prove statement (iv). First, let us write 

CQ{PoPi) = V{p.,e)nB{^i±^,^)nn^i'jHQ,K,.fe). 

Since Q is bounded, there exists po E Q such that 
CqiPijPj) C Q C B{po,diam{Q)). Therefore, we can write 

CQip^,Pj) = V(p„e) n ^) 

nrifei'i (^^Q.Ci-j^i.fc) nB(po,diam(Q))) . 

The map pi — » V(pi, e) is upper semicontinuous at any point 
Pi G Qe because there exists a neighborhood around any 
point Pi where the visibility set can only shrink. Also, the 
map Pi V{pi, e) has a compact range. Thus the notions of 
upper semicontinuity and closedness are equivalent and, hence. 
Pi V{pi, e) is closed at any point p,; G Q^. The second term 
is the closed ball of radius | centered at the point Ei^h. and 
clearly depends continuously on pi. 

Let oi, . . . , Ok be the strict concavities of Qe. Now Vij^k 
belongs to exactly one strict concavity, say a^^. Given a^^, 
Vi,j^k is a function of {pi,pj). We can write Wij.fc := 
Vij,k{pt,Pj) = argmin{dist(u, [pi,Pj]) \ v G a.^^}. Since a.,^ 
is a continuous curve, it is clear that Wij.fc is a continuous 
function of {pi,pj). We now show that HQ^{vij^k{Pi,Pj)) H 
B{po,diam{Q)) is in turn a set-valued map continuous with 
respect to {pi,pj). Since Vij^k is a continuous function of 
{pi,Pj), from Propositions 1 and 8 in [17], it suffices to show 
that Hq^{v) n B{po,diam{Q)) is continuous at any point 
V over the domain Ui^. This is equivalent to showing that 
^Qci"^) ^ B {po , di'Am{Q)) is closed as well as open at any 
point V G Qi^. Without loss of generality, let us assume that 
OHq^ (v) is parallel to the x axis as shown in Figure [17] 
Let us first show that Hq^{v) n i?(po, diam(Q)) is open at 
any point v G a^^. Let {vi} be a sequence in a^j. such that 
vi V, and let q be any point on dHq^ [v). Let qi be the point 
on dHQ^{vi) that is either vertically above or below q. Now 

k—'Zill — diam((5), where s{vi) is the slope of the line 

defining the boundary of Hq^{vi). Since Oi^. is continuously 
differentiable, as vi v, we have that s{vi) s{v) = 0. 
Thus, qi q. Hence, Hq^ n i?(po7 diam((5)) is open at any 
point V G ai^. We now show that Hq^ H B {po , dia.m{Q)) 
is closed at any point v G Ui^. As before, let {vi} be a 
sequence in a.jj, such that ^ v. Let G Hq^{vi) D 
B{pQ,dmm{Q)) and let q. But dist((7;, iJg^ (u) fl 

i?(po, diam(Q))) < \s{vi)\ dmm{Q) and s{vi) s{v) = as 



vi V. Therefore, dist{qi, Hq^{v) fl i?(poi diam(Q))) 0. 
Since qi ^ q and the distance of a point to the convex set 
-^Qs(^) ^ ^(Po: diam(Q)) is a continuous function, we de- 
duce dist(5z, iJQ^(u) n B(po, diam((3))) — > dist(g, i?Q, (u) n 
B(po,diam((5))) = 0. Hence q G i/g, (w)nB(po, diam(Q))). 
Therefore Hq^{v) fl i?(po: diam((3))) is a closed map and, 




Fig. 17. Illustration of various notions used in proving that Hq^ (v) n 
iJ(pOi diam((3)) is continuous at any point v over the domain ai^. The 
circle represents the boundary of the set S(poi diam(Q)). The arc represents 
ai^, . The lines tangent to a^^. at v and vi are the boundaries of the half- 
planes Hq^(v) and Hq^ivi) respectively. The interior of the half-planes 
is in the direction of the an'ows. In the figure on the right, & represents 
dist(g, , Hq^ {v) n B(po, diam{Q))). 

in turn, is continuous over the domain ai^. This implies that 

Hq^ {vi,jM{pi,Pj)) n i?(po, diam((3))) is continuous at pi for 
fixed Oij.. 

Now, let {pfjp'j) Pi and let q^ ^ q where q^ G 
CgiPiTPj)- We need to show that q G CQ{pi,pj). Since 
q'' G CQ{pf,p^), we have that q'' G V(pf,e). Since the map 
Pi V{pi,e) is closed, we have that q G V{pi,e). Also, 
g 5(£i±Pi^ r). Since the map {p„pj) -> B{S^, §) is 
continuous, it follows that q G ^ r). Therefore, what 

remains to be shown is that q G HQ^{vi j i{pi,pj)) for all I G 
{1, . . . , riy }. In what follows, for any points {pi,Pj), {Pi,Pj) 
and {p',,p'j), let Vij^i{pi,pj), Vi^j^i{p'y ,p'^) and Vi^j^i{p'i,p'j) 
belong to the concavities a^, , a*^ and a^^ respectively and with 
I belonging to {1, • • • {I,-- - ,ny} and {I,-- - ,n-j} 

respectively. Since (7*^ G HQ^{vi^j^i{p^ ,p^)), if there exists /cq 
such that for all k > fco, we have that a*j G {0,^ , • • • , a^^ . . }, 
then we are done. Equivalently, it suffices to show that if a is 
any strict concavity such that a ^ {ai-^ , a,i„ . . } then there 
exists a neighborhood B around pi such that for all p'^ G B, 
a ^ {a'^_^ , • ■ • , , }. But if a ^ {a^j , ■ • ■ , a^^ . , }, then a is 

in the exterior of at least one of V(pi,e), , |) and 

HQ^{vij,i{pi,pj)). The existence of S now follows from the 
fact that the maps pi — > V(pi,e), pi ^ jj|- p.+Pj ^ r-^ ^j^^ 
Pi — > HQ^{vij^i{pi)) Ci B{pQ,dia,mQ) are all closed. ■ 
Proof of Lemma \IV.5\ The second inclusion in fact (i) is 
a direct consequence of the definition of Q\c,g. We prove 
now that Semst.s Q Q\c,g by contradiction. Let V be any 
point set and assume, without loss of generality, that Q{V) 
is connected (otherwise, the same reasoning carries over for 
each connected component of Q{V)). For simplicity, further 
assume that the distances \\pk — Pi\\, k,l G {l,...,n}, k ^ I 
are all distinct. This ensures that there is a unique Euclidean 
Minimum Spanning Tree of GiV). If this is not the case, 
the same reasoning exposed here carries through for each 
Euclidean Minimum Spanning Tree associated with Q{V). 
Let {pi.pj) G feEMST.s(^) and ^ £g,,,giV). Since 
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necessarily {pi,pj) S £g{P), the latter implies that there 
exists a maximal clique V' of the edge {pi,pj) in Q such 
that {pi,Pj) ^ ^SemstI^')- If we remove the edge {pi,Pj) from 
GEMSxgi'P), the tree becomes disconnected into two connected 
components Ti and T2, with pi G Ti and pj G Now, 
there must exist an edge e G i^SemstI^') with one vertex in 
Ti and the other vertex in T2 and with length strictly less 
than \\pi — pj\\. To see this, let {ei, . . . , e^} be the edges of 
0emst(7'') obtained in incremental order by running Prim's 
algorithm (e.g., see [18]) starting from the vertex pi. Because 
Pi is in Ti and pj is in T2, there must exist at least an 
edge in {ei, . . . ,ed} with one vertex in Ti and the other 
vertex in T2. Let s G {1, . . . ,d} he such that Cs is the 
first edge having one vertex in Ti and another vertex in T2. 
Since {pi,Pj) ^ ^SemstI^')' according to Prim's algorithm, the 
length of Cs must be strictly less than \\pi —pj\\ (otherwise, the 
edge {pi , pj ) will be part of the Euclidean Minimum Spanning 
Tree of V'). If we add the edge e., to the set of edges of T1UT2, 
the resulting graph G is acyclic, connected and contains all 
the vertices V, i.e., G is a spanning tree. Moreover, since the 
length of Cs is strictly less than \\pi — Pj\\ and Ti and T2 
are induced subgraphs of 5EMST,e(^), we conclude that G 
has shorter length than GEMSxgiV), which is a contradiction. 
Fact (ii) is a consequence of fact (i). ■ 

Appendix II 
Additional analysis results and proof of 
Theorem IV.4I 

In this section we provide some key technical results that 
help establish Theorem IV.4I We first present the technical 
results on which the proof depends. 

Lemma II.l (Properties of the relative convex hull) For 

any allowable Q,Qe, the following statements hold: 

(i) if Vi,V2^ Q, then xco{Vi,Q) C rco(7'i U V2, Q); 

(ii) if rco{P',Q) is a strict subset of rco{P" ,Q), then 

(iii) ifV CQe and g{r) C g,^^,{V), then for p, G V, we 
have that {pi] n Vc(rco(7', Q^)) C Ve(co(M.g)). 

Proof: Statements (i) and (ii) are obvious from the 
definitions of the relative convex hull and its perimeter To 
prove statement (iii), assume pi ^ Ve(co(A/'i,g)). Then, pi 
belongs either to dco{Mi,g) \ Ve(co(A/i,e)) or to the interior 
of co{Mi,g). In the former case, pi G [pk,pi\ where Pk,pi G 
Ve(co(7\/i,£;)); see Figure [Tsl left). Note that since Q C C/sens 
and Pk, pi are neighbors of pi in Q, we have [pi,pk] Q Qc 
and [p^,pl] C Q^. Since p., G [pk.Pi], this implies that 
[PkiPi] C Qg. Therefore [pk^Pi] C rco{P,Q^) because the 
shortest path in between Pk^Pi is contained in rco(P, Q^)- 
Thus, Pi cannot be a vertex of tco{P, Q^). Let us now look 
at the case when pi belongs to the interior of co(A/i_g); see 
Figure [TSjright). Then, since the strict concavities of are all 
continuously differentiable curves, the shortest path between 
any two points Pk,Pi G Vc(co(A/i.g)) will be a sequence 
of straight lines and curves as shown in Figure [Tsl right). 
By the definition of the relative convex hull, the shortest 
paths between two points are contained in the relative convex 




Fig. 18. Illustration of the cases when pi Vc(co(A/i,g)). The point set 
represented by the black disks is Afi,g. The shaded region is co{Afi g). The 
straight line segments between any two points represent the fact that the two 
points are visible in Q^. The curved Hnes between pairs of points in the figure 
on the right represent shortest paths in Qe between the points. 

hull. Thus the region bounded by the shortest paths between 
adjacent vertices of co{J\fi.g) is contained in the relative 
convex hull. Clearly, since pi is the interior of this region, 
it cannot be a vertex of rco(7^, Q^). This concludes the proof 
of statement (iii). ■ 
See Figure [19] for a graphical explanation of Lemma III. 11 




Fig. 19. Illustration of Lemma Hi. II The outer polygonal environment is Q 
and the inner curved environment is Q^. The set of points represented by black 
and white disks are Vi and 7^2 respectively. Note that rco('Pi,QE), repre- 
sented by the dark shaded region, is a subset of rco('Pi \J'P2 , Qc), represented 
by the union of the dark and light shaded regions (c.f. Lemma lll. ll (i)). In fact 
rco('Pi , Qe) is a strict subset of rco('Pi U 7^2 , Qe) and hence rco CPi , Qe) 
has a strictly smaller perimeter than rco('Pi U7-'2, Qe) (c.f. Lemma lll. ll (ii')'). 
Now, Pi e Ve(rco(7'i WP2))- The set {pi} U {pi^, . . . ,pi^} is equal 
to the set Mi^g^^„^ where C?sens is such that r = +00 . Note that pi S 
Ve(co(A/i,gj^„j,)) where co{Mi,g^^^^) is the region bounded by the dashed 
lines (c.f. Lemma fll. 1 1 (iii)). 



Lemma II.2 (Properties of the constraint set) For any al- 
lowable Q,Q„ ifPeQ"^ and 02(7) C g^iP) C e.sens(P), 
then for all pi, the following statements hold: 

(i) The set Cp^,Q{J\fi,g.2)r\rco{Mi,g-i,V{pi,e)) is convex; and 

(ii) Cp,,Q(AAi,gi)nrco(M,ei, V(pi,e)) = Cp^, Q(7Vi, Si )n co(M,f; J. 

Proof Letp, g G Cp^, Q(A/i,e2) ri rco(A/i,ei , V(pi, e)). 
Since Cp.^Q{J\fi^g^) is an intersection of convex sets, it is 
convex as well. This implies that [p,q] Q Cp.^Q{J\fi^g2). 
Since Cp^^qlMi^g^) C Q^, we deduce that [p,q] C Q^. Now, 
p,q & rco(A/i,g;i , V(pi, e)). Therefore, by the definition of 
the relative convex hull, the shortest path between p and 
q in (5e is also contained in rco(A/i^£;i , V(pi, e)). I^^it the 
shortest path is the segment [p,q\ since [p,q\ C Q^. Thus, 
[p^q] ^ rco(A/i,ei , V(pi, e)). This completes the proof of 
statement (i). 

To prove statement (ii), note that rco(A/'i,gi , V(pi, e)) ^ 
co(A/'i,eJ. Hence Cp,,Q(M,ei) ^ rco{J\fi^g^,V{pi,e)) C 
Q'i-Q(-^j,ei)'~'<3o(-^j-ei)- To prove the other inclusion, let q G 
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Cpi,Q(M.ei)nco(M.gJ. Since q G co(7\/i,ei), then q belongs 
to the closed triangle formed by pi and two other points 
belonging to Yc{co{Ni.g^)), say Pk,Pi', see Figure [20l Also 
9 £ ^Pi,Qi^i,Gi)- From the construction of Cp-,Q{Mi^g^), it 
follows that q G V{pj, e) for all pj G Mi,g^. In particular, this 
means that C and C Qe. 

Since Pi,pi G M.ei and [pi.pi] C V(pi,e), we deduce that 
[PiiPi] ^ i'co(7\/i,g;i , ^(Pii £)) because the relative convex hull 
contains the shortest path in V(pi, e) between any two points 
contained in it. Let the line containing the segment [q^pk] 
intersect the segment [pi^pi] at q' . Note that q G [(?',Pfc] be- 
cause by choice of pk, pi, we have that q belongs to the closed 
triangle with vertices pi, pk, pi- Now, since q' G [pi,pi] and 
[pi,pi] C rco(M,ei, V(pi,e)), then q' G rco{N'i,g^,V{pt, e)). 
We now claim that the segment [q',pk] C V{pi,e). To see 
this, note that [q' ,pk] is a subset of the polygon given by the 
ordered set of vertices {pi,pk, q,pi)- Since all the edges of the 
polygon, [pi,pk], [Pk,q], [q,Pi], [pi,Pk], are a subset of Q^, 
we have that the interior of and the interior of the polygon 
{Pi,Pk-,q,Pi) do not intersect. Thus [q' ,Pk] C V{pi,e). Now, 
since q' ,pk G ico{Ni^g^,V{pi,e)), the shortest path between 
q' ,Pk in e) is also contained in rco(A/i^ej 7 ^{Pi^ c))- Thus 
q G [g'lPfc] ^ rco(A/i.gi , V(pi, e)). This completes the proof 
of statement (ii). ■ 




Fig. 20. The set of points represented by black disks is Mi- The shaded 
region represents co(A/i). q, represented by the white disk, is any point in 
co{J\fi). Since q is e-robustly visible from every point in Afi, the dashed line 
segments are a subset of Qe. Also every point in Afi \ {pi} is visible from 
Pi by the definition of Afi- Hence the solid line segments are also a subset 
of Q,. 

Statement (i) in the above lemma states that the constraint set 
is convex. Statement (ii) gives an alternate way of computing 
the set Xi in step 3: of the Perimeter Minimizing 
Algorithm when ^/constr and Q^ens are identical. 

In what follows, we analyze the Perimeter 
Minimizing Algorithm by means of the closed 
perimeter minimizing algorithm over any proximity graphs 
Q2 '^Gi Q Gaens, defined as follows: 



T, 



Ui , where 

P.,y(M,<?2) n rco(M,ei > V(p,;, e)), 
Pk^P, 



lim CC(X,(Pfc„J) 

kr,-,. —>--\-OC 



MP) 

Ki{P) = { 

{CC{X,{Pk^)} is a convergent subseq. of {CC(X,(Pfe)}| 

— (Pi - Pi) Pi ^ K'i{P) 



Ui G 



\\pi - Pi 



and where Tg g {P)i is the zth component of Tg g (P). 

Note that Tgl.g^X/') C T-g^^g^{P) if Gi = Ssens 
and G2 = ^/constr- Therefore, the evolution under the 
Perimeter Minimizing Algorithm given by the tra- 
jectory {P[i]}tez>o with P[t+l\ G Te^,„„e™.,(^^M) is just one 
of the possible evolutions under Tg g given by the trajectory 
{F^[t]}tez>o with P^[t+1] G T^lgl{P^[t\). We now begin 
stating the intermediate results for Tg_^ g^. The first of these 
results essentially states that for at least one robot located at a 
vertex of the relative convex hull, the set Xi is large enough 
for it to move under the application of Tg_^ g^. This, in turn, 
implies that the relative convex hull shrinks. 

Lemma II.3 (Nontrivial constraints) Let Q, Qf be allow- 
able environments. Assume that P G and the proximity 
graphs Gi and G2 satisfy 

A) not all points in P are equal; and 

B) 02 (P) C gi{P) C G.<,UP) and G2{P) is connected 
Then there exists pi G Ve(rco(P, Qe)) such that diam(Xi) > 
0, where X^ = Cp^^M,^^ H rco(M,gi , V(pi, e)). 

Proof: Note that pi C Xi and Xi is convex. Thus, 
diam(Xi) = if and only if Xi — {pi}. The remain- 
der of the proof is organized as follows. We first establish 
some necessary conditions for Xi = {pi}. We then show 
that under the assumptions in the lemma, those necessary 
conditions cannot hold for every pi G Ve(rco(P, Qe)). 
Note that X, = Cp.,Q(M,gJ n 1x0(^^5^ , V(p^, Qe)) and 
Cp.,Q(A/'i,e2) = n ^QiPi^Pj)- For Pj ^ -^^02^ 

Pj6A/'i,S2 

CqiPi.Pj) = V(K,e) n P( ^' ^ , ^) Pi HqXvi^jM) where 

fe=i 

Vi,j^k's are computed at step 3: of the Constraint Set 
Generator Algorithm. Thus, we can write 



nrco(M,Bj,V(p,,Q,)), 



(7) 



which can be rewritten as 



1 n ri ^^Q. I • 

Since rco(7\/i,ei , V(pi, e)) ^ V(Pije)> we have that 
V{pi, e) n rco(M,ei I ^{Pi^ Qe)) = rco(M,ei 1 e)). From 
assumptions (A) and (B), any robot i will have a neighbor 
in Gi{P) and G2{P) that is not placed at p.;. Hence, A/i.e^ 
and A/i,g;2 are point sets strictly larger than {pi}. Therefore, 
TCo{Afi^gi,V{pi,e)) is strictly larger than {pi}. Therefore, 
there exists a neighborhood around pi in rco(A/i.gi 1 ^(Kj 
as shown by the shaded region in Figure |23] Let this neigh- 
borhood be B. 



Therefore, if pi belongs to the interior of 

then Xi will be strictly larger than {pi}. Hence, if Xi = {pi} 
then Pi belongs to the boundary of one or more of the 
sets ^ ^) and HQ^{vij^k), where pj G M.Sa 

fc e {1, . . .,nij}. 

If \\pi ^ PjW < r, then pi belongs to the interior of 
^(^4^, §); see Figure |2T] (left). Also, if fe ^ then 
Pi belongs to the interior of HQ^{vi,j^k)\ see Figure |2T] (right). 
Therefore, if pi belongs to the boundary of HQ^{vi,j^k), 
then Vi j k G and if pi belongs to the boundary of 

r), then \\p, ~ p,\\ ^ r. 
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Fig. 23. The thick lines represent the boundary of rco(A/i,5j , V(pi, e)). 
The shaded region is a subset of rco{J\fi^g-^,V{pi,e)) and forms an entire 
neighborhood around pi. HQ^{vi j j) and HQ^{vi f^^) are the half-planes 
with boundary being the lines passing through [pi,Pj] and \pi,Pk\ respec- 
tively and interior in the direction of the arrows. In the figure on the right, 
WPj ~ Pill = ^ the angles between the vectors pj — pi and pk — Pi is 
greater than ^. The circle represents ^ 




Fig. 21. In the figure on the left, \\pi — pj\\ < R. The circle represents 
£j^Pi+P£ ^ ^). In the figure on the curved line represents a strict concavity. 
^i,j,k represents the point on the concavity nearest to [pi,pj\. It can be seen 
that dist{vi j k, \pi,pj\) > 0. The half-plane tangent to the strict concavity 
at Vi j k and with interior in the direction of the arrow is Hq^ {""i.j.k)- 

Since pi G Ve(rco(P, Q^)), from Lemma III. II (iii) 
we have that pi G Ve(co(A/i.g2)). This implies that 

B{ ^^ 2 ' '2} ^ convex set strictly containing 
{pi}; see Figure l22] 




Fig. 22. Pi is a vertex of the convex hull of the set {pi,Pj,Pk\- Here 
\\Pi ~Pj II = \\Pi ~Pk\\ = ^- Note that the angle between the vectors pj —pi 
and pk — Pi is strictly less than tt. The circles represent B{ ^^^^ , ^) and 
j^^ Pi+Pj ^ fc J ji^g shaded region represents their intersection. If \\pi — 
PjW < r or]|pi — Pfc II < r, the region of intersection will still strictly contain 

{P^}■ 

Therefore || B{-^—^ — ~'2^ ^ ^ strictly contains 

{pi}. Therefore, pi must belong to the boundary of some 
Hq^Vij^k)- It is easy to see that HQ^{vij^k) n S is strictly 
larger than {pi}. Therefore, pi must belong to the boundary 
of another set HQ^{vi k.i) or ^^ P'+Pj ^ ^s explained in 
Figure [23] 

Therefore, forpi G Ve(rco(P, Qt)), we obtain the following 
necessary condition for diam(Xi) ~ 0. There exist at least two 
neighbors of pi located at pj , pk such that one or more than 
one of following cases hold: 

(A) \\pj-pi\\ = r, Vi^kd e [Pi,Pk] for some ^ G {1, . . ■,nik}, 
the inner product {pj — Pi,Pk — Pi) < and pj ^ 

Hq^ivi^kiy, 



(B) G [Pi^Pj], Vi^k,i G [Pi.Pk] andpfc ^ iJQ,(vjj,m), 

Pj ^ HQ^(vi,k.i) for some m G {1, . . . and some 

/ G {1, . . . ,71ifc}. 

Cases (A) and (B) correspond to Figure |24] left and right 
respectively. 




Pi Pii 



Fig. 24. Illustration of cases (A) and (B) in the proof of Lemma lll.3l These 
cases illustrate the necessary conditions in order for the constraint set Xi to 
be a point. 

We now show that the necessary condition cannot hold 
for every pi G Ve(rco(P, Qe)). Let pi^ G Ve(rco(P, Qe))- 
Then pi^ has at least two neighbors in G2{P) located at 
Pj,Pk such that conditions (A) and/or (B) are satisfied. Now 
Pii G 9rco(P, Qe) because pi^ G Ve(rco(P, Qj)). Also, 
Vi^k,i G 9Qe n rco(P, Qe). Hence Vi^k,i G 9rco(P, Qj). Now, 
since does not contain any hole, rco(P, Q^) is simply 
connected. Thus the segment [pii,Vi,k.i] partitions rco(P, Qe) 
into two simply connected sets closed sets, and Q" such 
that n Ql ^ [p^, , v,,kj]- Let Pk & Q-,- ' 

Let us construct a segment by extending the 

segment [pi-^ , Vi^k.i] till it intersects the boundary of rco(P, Q^) 
at d. We refer to Figure |25] for an illustration of this argument. 
Then the segment [vi^k.i,d] partitions rco(P, Q^) into two 
components such that again there exists G Ve{P,Qe) in 
the component not containing pi-^. We point that = Vi^k,i, 
Pii ^ Pk and Vi^k,i = d are possible. Now, can only contain 
one neighbor in Q"^ and that can only be pi^ since all other 
points in Q"^ are not e-robust visible from . 

Let Pj be such that \\pj — Pi^W = r and the inner product 
{pj — PinPk — Pii) < 0. Let us construct a line segment [a, b] 
perpendicular to the segment [pii,Pj], passing through pj and 
intersecting the boundary of rco(P, Qg) at points a and b. 
The segment [a,b] partitions rco(P, Qe) into two components 
such that there exists p'/^ G Ve(P, Qe) in the component not 
containing p^^ . Here we point that p"^ = pj is possible. Now, 
p'/^ can only contain one neighbor in Q'^^ and that can only 
be pij since all other points in are at a distance strictly 
greater than r. 
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Fig. 25. Illustration of the construction of various segments and the partitions 
induced by them in the proof of Lemma lil.3l 

If, on the other hand, pj is such that there exists Wij,m G 
[pi,Pj] and pk ^ Hq^Vij^m), Pj i HQ^{vi^kd), then we pro- 
ceed in a way similar to the case of pk- Let us construct a Hne 
segment [vij^m,e] by extending the hne segment [pi,Wij,m] 
till it intersects the boundary of rco(P, Q^) at e. Then the 
segment [vij,m,s] partitions rco(P, Q^) into two components 
such that again there exists p"^ £ Vc(P, Qe) in the component 
not containing pi-^ . We point that p['^ = Vij^n, Pi\ = Pj and 
yi,j,m — e aie possible. Now, as before p"^ can only contain 
one neighbor in Q'^^ and that can only be j since all other 
points in Q[_^ are not e-robust visible from p"^. 

Thus, Xi^ = {pi^} implies that there exists a partition of 
rco(P, Qe) into Q-^ and Q'-^ containing p-^ € Ve(rco(P, Qe)) 
and S Ve(rco(P, Qe)) respectively such that p^_^ has only 
one neighbor in Q'/^ and vice versa. Now, let pi^ — p'/^. 
Therefore, pi^ ^ pi^. Again Xi^ = {^12} implies that 
rco(P, Qe) can be partitioned into Q'^^ and Q"^. Letp^^ G Q[^. 
Therefore, there exists p'/^ S Ve(rco(P, Qe)) such that pi^ 
cannot be a neighbor of p'l^ . Let pi, — p'l^ . Therefore, Xi^ = 
{pij} implies the existence of pi^ G Ve(rco(P, Qe)) such that 
Vii ^ {pinVi-i)- Let the cardinality of Ve(rco(P, Q^)) be n„. 
Proceeding recursively, Xi^ = {pi^ } implies the existence 
of Pi„^+i ^ {pii , ■ ■ ■ ,Pi„^ }■ This is a contradiction. ■ 
Now, let Gi(P) and G2(P) be graphs with fixed topologies 
and let G2 be a subgraph of Gi. Let ^01.02 = {P £ 

Q^|G2(P)CGi(P)Cg,e„s(P)}. 

We now present some results on smoothness of T^^ on 
^01.02 and contraction of Iperim.Q^ under the action of Tg^ g,^ . 
The technical approach in what follows is similar to the one 
in [7]. 

Lemma II.4 (Properties of the Perimeter 
Minimizing Algorithm) The set Ygi.G2 '^"'^ 
map Tq_^ have the following properties: 

(i) iGi.G2 compact, and Tq q is closed on Yg^.g^; 

(ii) for all i e {!,.!'. G2^P)i ^ 
(rco(P, Qe) \ Ve(rco(P, Q,))) U {pJ; 

(iii) ifG2 is connected, then there exists pi £ Ve(rco(P, Qe)) 
such that T^^_GSP)i ^ rco(P,Qe) \ Ve(rco(P, Qe)); 
and 

(iv) rco(Ta^,G^(P), Qe)) C rco(P, Qe) at any P e Yg,,g,, 
where Tq_^ {P)i is the ith component of Tq^ ^.^ (P). 

Proof: We begin by proving statement (i). Showing that 
YG1.G2 is compact is equivalent to showing that ^01,02 is 



closed and bounded. Since 1gi,G2 ^ then 101,02 is 
bounded. Now, let {Pk} be any sequence in ^01,02 such 
that Pfe P. Let us show that P e Ygi.g2- By definition, 
G2(P) C Gi(P). It remains to be shown' that Gi(P) C 
Gsens{P)- To see this, let {pi,Pj) be an edge that does not 
belong to Gsens{P)- We now show that it cannot belong to 
Gi(P). Then, either [p„pj] Q, or ||p, -p,|| > R. If 
[pi,Pj] <^ Qe, we can construct neighborhoods around pi and 
Pj such that for all p^ and p'j belonging to the neighborhoods, 
we have [p^,Pj] ^ Qe\ see Figure |26] If, on the other hand, 

® V V^ ® 

Pi ^ Pj 



Fig. 26. If [pi,Pj] (2 Qe, there exist neighborhoods of pi and such that 
[Pi'Pj] ^ Q«' ^" Pi Pj belonging to the neighborhoods. 

^ Pjll ^ R + S, then let p^,Pj be any points such that 

M^P.W < I and Hp;- -p,|| < |. Then U - p'^\\ = 

\\P^~P^+P^-PJ_ + P3-P'J\\ > \\P^ - pA ~ M - P^\\ - \\p'j ' 

PjW > P + (5 — I — I = R. Thus there exist neighborhoods 
around pi and pj such that for all p^ and p'j belonging to the 
neighborhoods, we have \\p^ ~ p'j\\ > R- Therefore, we can 
construct a neighborhood around P in Qe such that for all 
P' in the neighborhood, (p';,p^) is not an edge of ^sens(P') if 
{pi,Pj) is not an edge of ^?sens(P)- Since P^ P, there exists 
fco such that for all k > ko, Pk belongs to this neighborhood. 
Therefore, for all k > ko, {p\,p^) is not an edge of Qssns{Pk)- 
Since Gi(Pfe) C GsensiPk), for all k > ko, {p'^,Pj) is not an 
edge of Gi(Pfc). Then, since Gi(P) has the same topology as 
Gi(Pfc), we have that {pi,pj) is not an edge of Gi(P). This 
completes the proof that Ygi,G2 is compact. 

Showing that Tq_^ is closed over yG'i,G2 is equivalent to 
proving that the map P Ki{P) is closed over Ygi,G2 for 
each is {!,..., n\. Now, let Pi ^ P and let qi ^ q where 
qi e Ki{Pi). We need to show that q e /^^(P). Since qi € 
Ki{Pi), we have that qi = limm^oo CC(Xi(P^. )) for some 
sequence P^. Pi. Hence, given 5 > 0, there exists mo{l) 
such that for all m > mo{l), we have \\qi-GG{Xi{Pl^J)\\ < 
6. Also, since Pj. Pi, there exists ki such that for all k > ki, 
we have \\Pk — Pi\\ < S. Without loss of generality, let us 
assume that A:„jj,(;) > ki. Now, let P" = P^ and = 
CC(X,(P/')). Therefore, ||P/' - Pi\\ < 5 and 'jfi'' - qi\\ < 5. 
Since Pi ^ P and (/; ^ g, there exists such that for all 
I > lo, we have that \\Pi - P\\ < S and \\qi - q\\ < S. But 
IIP, - P|| = IIP, - P/' + P/' - P|| > IIP/' - P|| - IIP, - P", ||. 
Therefore, ||P/'-P|| < (5+ ||P, -P"/|| < 2(5. Hence P/' ^ P. 
Similarly, it follows that for I > lo- we have H?" — 9|| < "26. 
Thus, q. It then follows from the definition of Ki{P) 

that q G is:i(P). 

We now prove statement (ii). It suffices to prove that 
K,iP) C (rco(P,Qe)\Vc(rco(P,Qe))) U {pJ. Let q G 
K,{P). Therefore, q = lim^^^oo CC(X,(Pfc„)) where 
{Pkm } is some subsequence of Pk P. We begin by showing 
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that the map P ^ Xi{P) is closed at P. Note that Xi{P) = 
Cp.,Q(M,G2)nrco(M,Gi, V(p,,e)) and that Cp,,q(M.g2) = 
^PJ&^f^.G2'^Q{P^^P3)■ If G M,G2. then is an edge of 

G2{P) C ^/sensC^*)- This implies that {pi,Pj) G J- But from 
Lemma riV.2l (iv). the map {pi,pj) CQ{pi,pj) is closed 
over J. Therefore, the map P CQ{pi,pj) is closed over 
^Gi,G2 for sach pj G Mi,G2- Since G2 has fixed topology, 
we have that the set of indices {j \ pj G A/i.G2} constant 
over Ygj^g^. Thus, the map P Cp.^Q{J^i^G2) is closed 
over YGi.Gg since by Proposition 4 in [17], the intersection of 
closed maps is closed. 

We now show that rco(A/i.Gi , e)) = 
rco(7\/i,Gi, Qe)- Since V{pi,e) C Qe, the inclusion 
rco(7\/i,Gi, V(pi, e)) ^ rco(A/i,Gi , Qe) follows directly. Let 
qi, 92 G rco(M,Gi, V(pi,e)). Therefore, qi, q2 G V(K,e). 
By definition the shortest path in V{pi,e) between qi and 
q2 is contained in rco(7\/i^Gi i Vd?;, e)). Since Qe does 
not contain any holes, the shortest path between qi and 
q2 in V{pi,e) is also the shortest path between qi and 
92 in Qe- Therefore, by definition of a relative convex 
set, rco(7Vi,Gi : V(pi, e)) is convex relative to Q^. Since 
M,Gi ^ rco(M,Gi, V(pi,e)) and rco(7Vj,Gi , Q^) is the 
intersection of all relative convex set containing Mi.Gi, we 
have that rco(A/i.Gi , Qe) ^ i'Co(A/i,Gi , V(pt, e)). This proves 
that rco(M,Gi, V(K,e)) = rco(7Vj,Gi , Qe)- 

Now, since the shortest distances between two points in 
Qe change continuously as the position of the points, the 
relative convex hull of the points also changes continuously 
as a function of the points. Again the set { j \ pj G J\fi.Gi } is 
constant over yGi.G2- Thus, the map P i'co(A/i , Q^) 
is continuous over yGi.G2 hence, also closed. Finally, 
since Xi{P) is the intersection of closed maps, we have that 
P X,{P) is closed. 

Then from the definition of a closed map q G Xi{P). 
But X,{P) C rco(M,Gi,V(p„e)) = rco(M,Gi,Qe) 
and since A/i.Gi ^ {Pi I {l:---7^}}^ we have that 
rco(M,Gi,Qe) C rco(P,Q,). Thus, X,{P) C rco(P,Q,) 
and hence q G rco(P, Q^). Now, let q G Ve(rco(P, Qe)). 
We show that then diam{Xi{P)) = 0. Our strategy is as 
follows. We first show that if q E Ve(rco(P, Qe)), then 
dia.m{X i{Pk^)) — > 0. But that in turn implies that q = Pi 
and hence pi G Ve(rco(P, Qe)). But pi G Ve(rco(P, Qe)) 
and diam(Xi(P)) > together imply that diam(Xi(Pfc^ )) 
does not tend to zero. This is a contradiction. Therefore, 
diam(Xj(P)) = 0. 

We now start by showing that dia,m{X i{Pk^)) 0. The 
range of the map P Xi{P) is compact, the notions of upper 
semicontinuity and closedness are identical. Therefore, P 
Xi{P) is upper semicontinuous at any point P G Fgi,G2- It 
follows from the definition of upper semicontinuity that given 
any (5 > 0, there exists mo such that for all m > mo, we have 
X^iPkJ C U,ex,{p)Bix,S). 

Since Xi{P) C rco(P, Qe), g is a strictly convex point on 
the boundary of Xi{P). By a strictly convex point, we mean 
that there exists a tangent to dXi{P) at q which intersects 
dXi{P) at exactly one point. Let qi be the point on the 
boundary of L)x^Xi{p)B{x,S) such that — q\\ = S; see 
Figure |27] It follows that qi is a strictly convex point on 



the boundary of lJxeXi(P)B{x,S). Let li be the tangent at 
qi. Now, let I2 be the tangent to B{q,6) parallel to li. It is 
easy to see that as (5 0, the region of lJx^Xi{P)B{x,S) 
between the lines li and I2 tends to the point set {q}. Now, 




Fig. 27. Illustration of notions used in the proof of Lemma III. 41 Tlie shaded 
region represents Xi{P). The outermost boundary represents the boundary 
of Uj;gXi{P)^{^: The small circle represents B{q, 5). 

since GG{Xi{Pk^ )) q, there exists mi such that for all 
m > mi, CC{Xi{Pk^)) G B{q,S). For m > max{mo,mi}, 
let y be a point on the opposite side of I2 as CC{Xi{Pk^)) and 
equidistant to I2 as CC{Xi{Pk^)); see Figure [27] All points 
on the right of I2 are nearer to y than to CC{Xi{Pk^)). Also, 
the distance of y to any point in lJj.(=Xi{P)B{x, S) on or to the 
left of I2 tends to as 5 0. 

Let .Ti = argmax,{||2; - CC{X,{PkJ)\\ \ x G X^{PkJ} 
and let X2 = arg max:c { 1 1 a; - y|| | a; G Xi{Pk,^)}. By 
definition of the circumcenter and the fact that it is unique, 
\\xi — CC{Xi{Pk^))\\ < ||a;2 — y||. Also, by construction |[a;i — 
CC(X,(Pfe„))|| > \\x2 - CC(X,(Pfc,J)||. It follows then that 
|a;2~2/|| > ||x2 — CC(Xi(Pfc^))||. Therefore, a;2 belongs to the 
left of l2- But as pointed earlier, ||a;2 — — *■ as S —> 0. But 
the circumradius of Xi{Pk^J is equal to ||xi— CC(Xi(Pfe^))|| 
and \\xi —CC{Xi{Pk^))\\ < \\x2 — y\\- Since (5 can be chosen 
arbitrarily small, it follows that H^i - CC{Xi{Pk^))\\ 
as m 00. This proves that if q E Ve(rco(P, Qe)), then 
diam(X,(Pfe )) ^ 0. 

Now let be the ilh component of Pk^- Since 
CC{X,{PkJ) G X,{PkJ and diam(X,(Pr„J) ^ 0, it 
follows that WPi"- -CC{Xi{PkJ)\\ 0. But Pfc,„ ^ P, and 
therefore, p^"^ — > pi. This implies that CC{Xi{Pk^)) pi. 
Thus Pi = q and Pi G Ve(rco(P, Qe)). We now claim that if 
diam(Xi(P)) > and pi G Vc(rco(P, Qe)), then there exists 
a neighborhood of P in 101.02 such that for all points P' in 
the neighborhood, we have that diam(Xi(P')) > d for some 
d > Q. But this is a contradiction since we have shown that 

dia,m{Xi{Pk^)) — > for some sequence Pk^ P. 

To prove the above claim, note that since V{pi,e) C 
rco(M,Gi, V(pi,Qe)) and rco(M,Gi , V(pi, Qe)) = 
rco(M,Gi, V(pi,Qe)), we can write 

^. = rco(M,Gi,Q.)n I n B{^i±^.^)]n 

n ri ^Q.K.j,'.) I ■ w 

\PjeAf,,G2 '' = 1 / 

Now, let us first consider the case when Xi{P) has non- 
empty interior. Then one can choose a segment [p',p"] in the 
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interior of Xi (P) whose length is arbitrarily close to the length 
of diam(Xi(P)). Therefore, [p',p"] belongs in the interior of 
rco(M,Gi,Q£), B{^i^, §) and Hq^Vi^j^k) for all i, j and 
k. As shown earlier in the proof of statement (ii), the map 
Pi rco(7\/i,Gi (P), Qe) is continuous for every i. Therefore, 
there exists a neighborhood Bi of Yq^.g-i around P such that 
for all P' G Bi, the segment [p' ,p"] C rco(A/;;,Gi (P')> QO- 

From the proof of Proposition IIV.2I (iv), we know that 
the map {pi,pj) is continuous over J. It 

follows then that for fixed i and j, there exists a neigh- 
borhood of J around {p.i,pj) such that for all {p^,p'j) in 

the neighborhood, [p',p"] C B{ ^'t^' , |). Since G2 is fixed 
the set {j I 7 G A/i,G2} is fixed. Thus, there exists a 
neighborhood, S2 of Ygi,G2 of P such that for all P' G S2, 
we have [p',p"] C (np,eAr.,o, ^(^, §))■ Also, from the 
proof of Proposition IIV.2I (iv), for fixed j and a fixed strict 
concavity a, the map {pi,pj) — > HQ^{v{pi,pj)) is continuous 
over J, where v{pi,pj) is the point on the strict concavity 
a nearest to the segment [pi,pj]. Again, from the proof of 
Proposition IIV.2l (iv). there exists a neighborhood of J around 
(pi , pj ) such that for all points {p'^ , p'j ) in the neighborhood, we 
have a'^^ G {a^^ , • • • , a^^ . . } for any Z G { 1 , . . . , n'^j } where the 
notation is as in the proof of Proposition IIV.2I (iv). Let a'^^ = 
Ui^. Since belongs to the intenor of HQ^{vij,s{pi,Pj)), 

which in turn varies continuously as {pi,pj) for fixed a^^, 
there exists a neighborhood around {pi,Pj) such that for 
all (K,p;), we have [p',p"] C Hq^{v,,j^i){p'^,p'^). It then 
follows that there exists a neighborhood B3 of Ygi.G2 
around P such that for all P' = {p[, ■ ■ ■ ,p'j^) G B3, 

[p',p"] c (n,^,^^^^(p,)nr!fii^Q.K.,fc(K,p;)))- Thus, 

for any P' G 61 nSz n^s, we have that [p',p"] C Xi{P') or 
that diam(Xi(P')) > d where d is equal to the length of the 
segment [p',p"]. 

Second, let Xj (P) have empty interior. Therefore, Xi (P) is 
a line segment and since pi G Ve(rco(P, Qe)) we have that pi 
is one of the end points of Xi{P). Also, since Xi{P) is a line 
segment, it follows from equation (|9]l that Xi{P) must belong 
to the intersection of the boundaries of at least two of the 
sets rco(7\/i,Gi , Qd) '^he half-planes -ffg^ (wij.fc). The set 
Xi{P) cannot belong to the boundary of a ball P( ^'"^^^ , |) 
since its boundary is curved. Therefore, Xi (P) belongs to the 
intersection of the boundaries of either (i) rco(A/i_Gj , Qe) and 
a half-plane HQ^{vij^k) or (ii) two half-planes HQ^{vij^k) 
and Hq^ivi^i^m)- Let us assume without loss of generality 
that Xi{P) belongs to the interior of the remaining sets. Then, 
as shown earlier, there exists a neighborhood, B4,, of Ig^ 
around P such that for all P' in the neighborhood, Xi{P) 
belongs to the interior of the remaining sets. Now cases (i) 
and (ii) correspond to Figure |28] top and bottom respectively. 
It follows immediately from the figure that we can choose a 
neighborhood, B5, of P in yGi.G2 such that for P' in the 
neighborhood, we have that the rco(A/i.Gi i Q<d ^ ^Qe {f^i.j,k) 
and HQ^{vi j k) H HQ^{vij^rn) contain a segment of length 
arbitrarily close to the length of Xi{P). Therefore, for all 
P' G S4 n 65, we have that Xi{P') contains a segment of 
length d where d > 0. 




Fig. 28. Illustration of cases (i) and (ii) in the proof of Lemma [11.41 The 
curved arcs represent strict concavities. The solid line segment represents 
Xi{P). The dashed lines represent the boundary of rco(P, Qe). Note that 
there must exist robots pr, Pa S ■^i,Gi ^™ the dashed lines to be the boundary 
of rco(P, Qe). Note that pi S Yc{RC H {P, Qe)) and also pi is an end point 
of the segment Xi(P). In the top figure Vi j jf , denoted by the white disc, 
is the point on the strict concavity a^j, neai'est to the segment [pi,Pj]. Here 
Pj G J^i,G-2- '^^^ half-plane Hq^{vi jj^) has interior in the direction of 
the arrow. In the bottom figure, Vi i „i and fi.j.fc are points on a^j and a^^. 
nearest to the segments \pi,p{\ and [pi,pj\ respectively. The direction of 
the arrows points towards the interior of the half-planes HQ^{vi i „^) and 
Hq, (I'ij.fc). Here pj, pi G M,G2- 

This completes the proof of the fact that if g G 
Vc(rco(P, Qe)), then diam(Xi(P)) = 0. But then this impUes 
that q ~ Pi. Thus, this completes the proof of statement (ii). 

To prove statement (iii), note that in proving state- 
ment (ii) we have shown that if diam(Xi(P)) > 0, then 
q = limfc,„_ooCC(X,(Pfc„J) ^ Ve(rco(P,Q,)). From 
Lemma HOI there exists at least one pi G Ve(rco(P, Qe)) such 
that diam(Xi(P)) > 0. Statement (iii) now follows directly. 

We now prove statement (iv). It follows from State- 
ment (ii), that PgiG2(-^) - rco(P, Qe). This implies that 
rco(Ta^,G^(P), Qe) C rco(P, Q,). U 

From Lemma III.4t i). we know that, in particular, if the 
graphs Qsen^i , Gmnsa havc fixed topology, then the closed 
perimeter minimizing algorithm is a closed map. However, 
during the evolution of the system we expect the graphs 
Gsens, Sconstr to switch. To study the convergence properties of 
"^Gsa^-Gcomu over such a dynamic topology, we define a set- 
valued algorithm which essentially embeds all possible evolu- 
tions that can happen due to switching. Given any allowable 
Q, Qe, define the perimeter minimizing algorithm over all 
allowable graphs to be the set-valued function T that maps 
points in Q" to all possible subsets of Qe by: 

r(p) = {P' G t^g,,gAp) I G2 C Gi C ese„s(P), 
Gi , G2 have fixed topologies 
and same connected components as t/sens}- 

Also, let 3^ = {P G Q" I 5sens(P) is connected} denote the 
set of all configurations where t/sens is connected. 

Proposition II.5 (Perimeter minimizing algorithm over all 
allowable graphs) The set y and set-valued map T have the 
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following properties: 

(i) y is compact, T[y) C y, and T is closed on y; 

(ii) rco(r(F), Qe) C rco(P, Q^), for P € Q^; and 

(iii) if P' e T(P) fl«t/ P e g^, then Fperim,Q.(^') < 

Proof: Since 3^ C Q", we have that y is bounded. Now, 
let {Pk} be a sequence in y such that Pk P. Then, from the 
proof of Lemma llL4l (i). we know that there exists fco such that 
for all k > fco, {p^,p^) is an edge of Qsens{Pk) implies that 
{pi.Pj) is an edge of f/sens(^')- But C/sens(^'fc) is connected. 
Therefore, Qsen^{P) is also connected. Hence, P € y and, 
therefore, y is closed. This proves that y is compact. 

Now let P' e T{P) where P ey. Then P' £ T^,,g„{P) 
for some Gi and G2 with fixed topologies having the same 
connected components as Q^snaiP)- Since P G 3^, 0sens(-P) is 
connected and hence, so are Gi, G2. The fact that P' 6 3^ or 
that Gsen>i{P') has one connected component can be verified 
exactly along the lines of the proof of Theorem |V.4| (i). This 
proves that Tjy) cy. _ 

Now, let {Pfc} be any sequence with P^ G T{Pk) such that 
Pfe — > P. Since Pk — > P, and the number of distinct graphs 
possible over n nodes is finite, there exists a subsequence 
Pfe„ ^ P such that Pfc„ G T^,^G,iPkJ where Gi(Pfc„) 
and G2(Pfe„) have fixed topologies for all m. Now, since 
-Pfc™ ^ -P' ^fc™ ^ ^ and Pgi^s closed over Yci.Ga from 
Lemma Iir4l (i), we have that P G Tgi.g,!^)- Now, to finish 
the proof, we have to show that G2(P) C Gi(P) C &ens(P)- 
The first inclusion follows from the fact the by definition of 
T(P) where G2(P) C Gi(P). Now, notice that we have 
shown in proof of Lemma III. 41 (i) that there exists fco such 
that for all fc > fco, {Pi,Pj) being an edge of GsensiPk) 
implies that {pi,Pj) is an edge of Gsens{P)- Also, by definition, 
Gi{Pkm) ^ 5sens(-ffc„)- Therefore, for all fc^ > fco. we have 
that (p^"'" , p^" ) being an edge of Gi {Pk,„ ) implies that [pi , pj ) 
is an edge of ^sens(P')- Therefore, Gi(P) C ^/sens(P')- This 
completes the proof of statement (i). 

Statement (ii) follows directly from Lemma III. 41 (ii). Fi- 
nally, statement (iii) is a consequence of statement (ii) and 

Lemma |in](ii)- ■ 
We are now ready to present the proof of the main result in 
the paper The proof uses the analysis results presented in this 
section and the discrete time LaSalle Invariance Principle for 
set- valued maps [7]. 

Proof of Theorem \V.4\ 
We begin by proving statement (i). Let [^o] , [^o] be 
two robots belonging to the same connected component of 
&ens(-P[io])- Then from Corollary IIV.6r ii). they must be- 
long to the same connected component of Gcon^a{P[to\)- 
Therefore, there exists a path between pi[^o] and Pj[to\ in 
GconA,{P[to])- Let {pk[to\,Pi[to\) be any edge of 0constr(-P[io]) 
that belongs to this path. Note that Uk[ti^] and w;[io] in 
the Perimeter Minimizing Algorithm algorithm are 
constrained to belong to the sets Cpf^^Q{Mj^g^^) - Pfe[io] 
and Cp^.Q{JVj.g^^^^J -Plito], respectively. Therefore, pk[to + 
l],Pi[to + 1] G CQipkhlPiM C n i?(MM+£iM, r). 
Hence, the edge {pj[tQ + l],pk[to + 1]) is an edge of 
Gsens{P[to + !])• Thus, there exists a path between Pi[to + 1] 



and Pj[to + 1] in Gsens{P[to + 1])' which proves statement (i). 

For statements (ii) and (iii), it suffices to prove the results for 
the system evolving under the action of T, since the evolution 
under Tg^,^^ is just one of the possible evolutions under 
T. Statement (ii), therefore, follows from Proposition III. 5f iii). 

We now prove statement (iii). From statement (i), it follows 
that the number of connected components of Gsens{P[t]) is 
nondecreasing. Since the number of possible connected com- 
ponents of Gsens is finite, this implies that, after a finite time t, 
the number of connected components of Gsens does not change 
and robots belonging to different connected components never 
detect each other's presence. To prove statement (iii), it now 
suffices to show that any two robots belonging to the same 
connected component of G^aK converge to the same point. 
Since the evolution of any group of connected robots from 
time t on is independent from the existence of other connected 
components, we can then assume, without loss of generality, 
that tJsens(Pp]) has just one connected component. 

According to Lemma [V3] and Proposition Ill.Sr iii). V^erim.Q^ 
is continuous and non-increasing along T on Q". From 
Proposition Ill.Sr i), T is a closed map on the compact set 
y C Q". Then, by the LaSalle Invariance Principle for closed 
set-valued maps [7], P[t] M, where AA is the largest 
weakly positively invariant set contained in 

{P G y I 3P' G r(P) s.t. Vi™.Q.(-P') = Vi™.Q.(^)}- 

Let us define the set diag(Q") = {(p, ...,p) G Q" | p G 
Qc} C y, and show that M = diag((5"). Clearly, 
diag(Q") C M.. To prove the other inclusion, we rea- 
son by contradiction. Let P G A4 \ diag((5"). Then, 
Lemma III.4I (iii) implies that for any connected G2(P), 
Gi(P) satisfying G2(P) C Gi(P) C Ssensl-P), there exists 
p,; G Ve(rco(P,Qe)) such that T^^ g^{P\ C rco(P,Q,) \ 
Ve(rco(P, Qe)). Therefore, if P' G T(P), then P' contains 
a number of robots N{P') belonging to Ve(rco(P, Qj)) that 
is strictly smaller than N{P). In turn, this implies that after 
at most N{P) steps, at least one vertex of rco(P, Qe) will 
contain any robot. By the definition of a vertex of rco(P, Qc), 
it follows that the relative convex hull of the configuration 
after at most N{P) steps is a strict subset of rco(P, Q^). 
From Lemma lll.lf ii) in Appendix [III V^nm,Q^ has strictly 
decreased. This is a contradiction with the fact that M is 
weakly positively invariant. 

To finish the proof, we need to show that every trajectory 
converges to an individual point, i.e., P[t] P* G Ai 
for all initial conditions pfo] G y. Note that {P[t]} is 
a sequence in a compact set. Hence, it has a convergent 
subsequence {P[tfc]}, say to P*. Since P[t] M, it 
follows that P* G A^. Using again Lemma IV. 31 we know 
that Vjjerim.Qs is continuous and Vpaim.Qc (P*) = 0. Hence, 
Vperim.Q,(-P[ifc]) Vpa-im.QAP*) = 0- This implies that, 
given any S > 0, there exists fco such that for all fc > fco, 

IVi.n.Q.(^[^fe]) - Vi.n,Q. ^ Vi.n,Q.(^M) < 

From Proposition Ill.Sl iii), we know that VJ,erim,Q^ (^'[i]) ^ 
^erim.Q, (-P[ifco]) for all t > tfcg . In turn this implies that 

'a set M is weakly positively invariant with respect to a map T onW if, 
for any uiq G W, there exists w € T(too) such that w £ W. 
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^perim.Q.(PM) < Vim,Q.(^[ifco]) < ^ for all <> tfc„ . HoW- 

ever, \\P[t] - P*\\ < nmax{||p,[i] - p*\\ \ i G {l,...,n}} < 
'^Vj,erim,Q, (^'[i]), wherc the second inequality holds because 
p* e rco(P[t], Qe), for i G {1, . . . , n}, and because the length 
of the shortest measurable curve enclosing a set of points 
is greater than the distance between any two of the points. 
Thus, for all t > tk,,, we have \\P[t] - P*\\ < nd. Hence 
P[t] P* eM. U 

Appendix III 
Proofs of results in SectionIvII 

Proof of Proposition \VI.1\ We begin by proving statement (i). 
Note that step 1: of the Constraint Set Generator 
Algorithm is in fact the sensing operation that takes 0{M) 
time and the result is a star-shaped polygon with at most M 
vertices. It is easy to see that the intersection of a star-shaped 
polygon with a half-plane takes time linear in the number of 
vertices of the polygon. Therefore, step 4 : of the algorithm 
takes at most AI time. If Q has at most k concavities, 
then the computational complexity of Constraint Set 
Generator Algorithm is 0{M)+kM, which is 0{kM). 

We now prove statement (ii). Step 1 : of the Perimeter 
Minimizing Algorithm is the result of the sensing op- 
eration taking 0{M) time. Since the sensor takes at most 
M readings, we can assume that the cardinality of A/i,e,„s 
is at most M. In step 2 : , Afi^g^^ is the set of all robot 
positions obtained in step 1 : . Let the computation of A/i,0„„,„ 
given the set Afi.g^„^ take r(M) time. Here t{M) depends 
on the specific heuristic used to compute the maximal cliques 
of a graph. In step 3:, Cp;,Q(A/i,eco„,i,) ^ '■^^ 

intersection of at most M convex constraint sets. It is clear 
from the proof of statement (i) that each of the convex sets 
CQ{pi,pj) has at most M vertices and can be computed in 
0{kM) time. The intersection of AI convex polygons of AI 
vertices takes 0(M^ logM) time [19]. Thus, the computation 
of Cp.,g(M,e„„,J takes MO(kM) + 0{AP logM) time and 
^Pi,Q{-^i,Gcoma) has at most vertices. The computation of 
rco(A/'i,as=n.> V(p^,e)) takes 0(M log M) time (cf. [20]), since 
the cardinality of both Afi^g^^^ and Vc{V{pi,e)) is at most 
M. Now, finally step 3 : involves the intersection of two star- 
shaped polygons Cpi,Q(A/i,e,o„sir) rco(7\/j,g,„, , V(pi, e)) 
having at most and M vertices respectively. Therefore 
the total number of vertices of the two polygons is at most 
+ M. The computational complexity of computing the 
intersection of the two star-shaped polygons is 0((Af^ + 
M) \og{AP + M) + M^) or 0{M^), where AP is due to the 
fact that the maximum number of possible intersection points 
is [21]. Thus the computational complexity of step 3 : is 
MO{kM) + 0{AP logA/) + 0(M) + 0((M2 + M) log(Af2 + 
M) + M^). Since k g 0{M), we have that the computational 
complexity of step 3 : is 0{M^). In step 4 :, the computation 
of the circumcenter of Xi, which is a convex polygon of 
at most AP vertices, takes 0{Al'^ log AI^) or 0(M^ log M) 
time [19]. Step 5 : can be performed in constant time. Thus, 
the complexity of the entire algorithm is 0{AI) + t{M) + 
0{AP) + 0{AP log a/) or t(A/) + 0{AP log A/). 



To prove statement (iii), note that if Q = ^/constr = ftens, 
then t{AI) E 0{1). Also, from Lemma III. 21 (ii), we have 
that Xi ~ Cp-^QlMi.g) n co(A/i,g). The convex hull of 
A/ points can be computed in 0{M \ogAd) time (cf. [19]) 
and hence co{Ni^g) can be computed in 0(A/logAf) time. 
The set Xi is an intersection of convex sets and can be 
computed in 0{AP + AI) or 0{AI^) time, and contains at 
most 0(Af2) vertices. Step 4: now takes 0{AP\ogM'^) 
or O(Af^logAf) time. Thus, the complexity of the entire 
algorithm is 0{AP log Af ). ■ 



